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Four  sets  of  boundary  value  problems  for  the  equations  of  fluid 
motion  and  transport  diffusion  are  solved  to  obtain  expressions  for  the 
spatial  variations  of  velocity  and  concentration  in  free  surface  flows. 
The  four  problems  investigated  treat  unidirectional  and  oscillatory 
uniform  shear  flows  in  rectangular  channels  of  infinite  and  finite  width. 
The  solutions  are  then  used  to  determine  longitudinal  dispersive  mass 
transport.  Analytical  expressions  for  the  longitudinal  dispersion 
coefficient  associated  with  infinitely  wide  channels  are  presented  as 
functions  of  the  temporal  and  spatial  maximum  of  velocity,  channel  depth, 
period  of  oscillation  (oscillatory  flow  only),  and  depth  mean  vertical 
coefficient  of  eddy  diffusivity.  For  channels  of  finite  width,  analytical 
expressions  for  the  longitudinal  dispersion  coefficient  are  presented  as 
functions  of  the  temporal  and  spatial  maximum  of  velocity,  channel  depth, 
channel  width,  period  of  oscillation  (oscillatory  flow  only),  and  cross- 
sectional  mean  values  of  vertical  and  lateral  coefficients  of  eddy  dif- 
fusivity. 

It  is  found  that  for  infinitely  wide  channels  the  longitudinal  dis- 
persion coefficient  in  oscillatory  flow  is  described  by  a  type  of  resonant 
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interaction  between  the  period  of  oscillation  and  the  time  scale  of  vertical 
mixing.  The  maximum  longitudinal  dispersive  mass  transport  occurs  when 
the  ratio  of  the  vertical  mixing  time  to  the  period  of  oscillation  equals 
1.58.  This  functional  behavior  is  shown  to  be  unique  and  distinctly 
different  from  the  dispersive  process  in  unidirectional  flow. 

For  rectangular  channels  of  finite  width  it  is  shown  that  the  dis- 
persive mass  transport  process  is  symmetric,  such  that  a  channel  geometry 
skewed  in  width  produces  the  same  dispersive  mass  transport  as  a  channel 
geometry  skewed  in  depth,  provided  that  specific  requirements  regarding 
the  maximum  velocity,  period  of  oscillation,  and  vertical  and  lateral 
mixing  times  are  satisfied. 

The  longitudinal  dispersion  coefficient  for  oscillatory  flow  in 
rectangular  channels  of  finite  width  is  presented  as  a  family  of  resonance 
curves  described  by  the  resonant  interaction  of  the  period  of  oscillation 
with  the  vertical  and^  lateral  mixing  times. 

For  unidirectional  flow  in  channels  of  finite  width  it  is  found 
that  the  predicted  non-dimensional  longitudinal  dispersion  coefficient 
increases  without  bound  as  the  skewness  of  the  channel  geometry  is  corres- 
pondingly increased.  Thus,  for  unidirectional  flow  the  three  dimensional 
solution  does  not  collapse  in  the  limit  to  the  two  dimensional  case  as  is 
found  for  oscillatory  flow. 

Predicted  values  for  the  longitudinal  dispersion  coefficient  are 
compared  with  values  based  upon  field  and  laboratory  data.  The  compari- 
son shows  that  when  reliable  measurements  are  available  for  all  of  the 
functional  parameters  included  in  the  predictive  solutions  the  agreement 
between  observed  and  predicted  values  is  good. 
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CHAPTER  I 
INTRODUCTION 

I-A  General  Background 

A  growing  awareness  in  recent  years  of  the  need  for  a  rational  approach 
to  the  planned  use  and  development  of  coastal  regions  has  significantly 
contributed  to  a  concern  for  maintaining  an  acceptable  level  of  water 
quality  within  our  river,  estuarine,  and  bay  system  waters.  To  achieve 
this  a  fundamentally  sound  understanding  of  the  physical  processes  which 
govern  the  mass  transport  characteristics  of  these  waters  is  essential. 

The  mass  transport  of  a  substance  as  it  applies  here  may  be  defined 
as  "the  time  averaged  net  transfer  of  substance  along  the  principal  flow 
axis  of  the  watercourse  or  channel."  With  the  aid  of  this  definition  one 
can  intuitively  reason  that  the  mass  transport  mechanism  for  a  conservative 
substance  is  governed  by  the  kinematic  description  of  the  flow  field  and 
the  distribution  of  the  substance  within  the  flow  field.  Unfortunately, 
the  complexities  associated  with  turbulent  shear  flows  in  natural  water- 
courses have  thus  far  precluded  a  complete  analytical  description  of  either 
the  flow  kinematics  or  the  distribution  of  the  substance.  However,  the 
analytical  treatment  of  longitudinal  mass  transport  in  steady  uniform 
axisymmetric  shear  flows  by  Taylor  (1,2)  provided  for  the  first  time  a 
reasonable  method  for  predicting  these  phenomena. 

In  his  analysis  Taylor  made  the  distinction  between  three  mechanisms 
of  longitudinal  mass  transport: 

(1)  Convective  transport  of  the  cross-sectional  mean  concentration 

of  substance  by  the  mean  cross-sectional  velocity, 
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(2)  Convective  transport  due  to  the  correlation  over  the  cross 
section  of  the  spatial  variations  from  the  mean  of  velocity  and 
concentration, 

(3)  Diffusion  of  the  substance  along  the  principal  flow  axis  (mole- 
cular diffusion  in  laminar  flow  and  eddy  diffusion  in  turbulent  flow). 

The  first  of  these  mechanisms  has  long  been  recognized  as  a  mass  transport 
phenomenon  and  is  the  most  amenable  of  the  three  to  an  analytical  treat- 
ment. In  the  absence  of  diffusive  transport  due  either  to  molecular  or 
turbulent  motion  this  mechanism  may  be  visualized  as  the  transport  along 
the  principal  flow  axis  of  a  slug  of  substance  uniformly  mixed  over  the 
flow  cross  section  and  undistorted  by  boundary  shear.  The  third  mechanism, 
that  due  to  longitudinal  diffusive  transport,  was  shown  by  Taylor  to  be 
less  than  1  per  cent  of  the  second  which  he  defined  as  dispersive  mass 
transport.  It  is  this  second  mechanism,  longitudinal  dispersive  mass 
transport,  which  is  the  subject  of  the  present  work. 

Following  Taylor's  initial  formulation  of  the  dispersion  mechanism 
investigators  began  to  extend  his  analysis  to  unidirectional  two  dimen- 
sional open  channel  flows  for  infinitely  wide  channels  and  subsequently 
to  the  more  realistic  three  dimensional  case  of  unidirectional  flow  in  a 
channel  of  finite  width.  More  recently,  the  emphasis  of  the  research 
effort  has  shifted  from  steady  unidirectional  flow  to  oscillatory  flow 
in  an  effort  to  apply  the  knowledge  obtained  from  the  former  to  develop 
a  predictive  capability  for  dispersive  mass  transport  for  the  latter.  A 
brief  summary  of  previous  efforts  follows. 


I-B  Discussion  of  Previous  Work  on 
Dispersive  Mass  Transport 

I-B.l  Unidirectional  Flow 

As  mentioned  previously,  Taylor  (1)  was  the  first  investigator  to 
define  the  mechanism  of  dispersive  mass  transport.  In  his  paper  Taylor 
considered  the  case  of  a  passive,  conservative  substance  transported  by 
a  laminar  steady  flow  in  a  circular  tube.  A  passive  substance  is  one 
which  is  present  in  small  enough  concentrations  so  as  not  to  alter  either 
the  velocity  flow  field  or  the  magnitude  of  the  diffusion  coefficient;  a 
conservative  substance  is  one  which  has  a  fixed  total  mass  present  for  all 
time  in  the  system  being  considered.  All  substances  referred  to  in  the 
present  work  should  be  considered  to  be  both  passive  and  conservative. 
To  begin  his  analysis,  Taylor  transformed  the  Eulerian  form  of  the  trans- 
port diffusion  equation  to  one  with  the  origin  moving  with  the  cross- 
sectional  mean  velocity  by  defining 

^  =  X  -  ut  (I-l) 

u  =  G  +  u"  (1-2) 

where  u  is  the  cross-sectional  mean  velocity,  u"  is  the  variation  from 
the  mean  of  the  velocity  in  the  radial  direction  due  to  the  shear  profile, 
and  X  is  the  coordinate  along  the  axis  of  the  tube  positive  in  the  direction 
of  flow.  Neglecting  molecular  diffusion  along  the  axis  of  the  tube  he  was 
then  able  to  write 

fl£^  +  ^<<  ^  =  JL  J-   fp^>  (1.2) 


where. 


3c 
(■^)>.   =  time  rate  of  change  of  concentration  as  seen  from  moving 

ot  c, 

coordinate  system 
=lc  +  -  9c 

at  "^  3x 


D  =  molecular  diffusion  coefficient  in  radial  direction, 
r 

If  (1-3)  is  integrated  over  the  cross-sectional  area  of  the  tube 
and  the  concentration,  c,  decomposed  into  its  cross-sectional  mean  and 
variation  component,  by  the  expression 

c  =  c  +  c"  (1-4) 

then  (1-3)  reduces  to 

(|f)5  =  -^I^^  (1-5) 

For  standardization  throughout  the  text  overbars  shall  be  considered 
to  represent  cross-sectional  averaging  whereas  the  symbol  <  >y  shall  be 
considered  to  represent  temporal  averaging  over  the  periodic  interval  T. 
The  form  of  (1-5)  led  Taylor  to  define  the  convective  transport  shown  as 
a  Fickian  flux  which  he  called  dispersion,  or 

E,  =  -  -!—  \\   u"c"dA  (1-6) 

A^  ^^A 

where  E,  is  the  longitudinal  dispersion  coefficient  for  steady  unidirec- 
tional flow  and  A  is  the  flow  cross-sectional  area.  As  indicated  by  (1-6) 
the  heart  of  Taylor's  analysis  lies  in  the  determination  of  the  spatial 
variations  of  velocity  and  concentration,  u"  and  c"  respectively.  The 
velocity  variation  was  readily  determined  from  the  parabolic  velocity 
profile.  To  determine  c"  Taylor  returned  to  (1-3),  decomposed  c  by  (1-4) 
and  made  the  following  assumptions: 

(1)  Steady  state  conditions  exist  with  respect  to  the  moving  coor- 

3c 
dinate  system,  i.e.  (g^V  ~  ^• 

(2)  c"  is  not  a  function  of  ^. 

3c 

(3)  ^F  is  constant, 

dt, 


These  assumptions  reduced  (1-3)  to 


D 

..  AC 

u 


85    r  3r  ^  9r  ''  ^       ' 

which  Taylor  then  solved  for  c",  performed  the  integration  of  u"c"  in 
(1-6)  and  produced  the  following  expression  for  the  longitudinal  disper- 
sion coefficient 

a2u2 


where  a  is  the  radius  of  the  tube,  and  u„^„  is  the  centerline  velocity. 

maX 

It  is  important  to  note  that  the  dispersion  coefficient  varies  directly 
with  the  ratio  a^/D  which  represents  the  time  scale  of  mixing  of  the 
substance  over  the  flow  cross  section  due  to  transverse  diffusivity. 
Using  the  definition  given  in  (1-6)  for  E^,  (1-5)  may  then  be  written  as 


Equation  (1-9),  with  E,  given  by  (1-8),  shows  that  the  cross-sectional  mean 
concentration  of  a  substance  is  effectively  diffused  longitudinally  about 
a  point  traveling  with  the  mean  velocity  by  the  combined  effect  of  the 
velocity  shear  profile  and  the  transverse  diffusivity.  Taylor  correctly 
cautioned,  however,  that  in  applying  (1-9)  the  time  scale  of  convective 
transport  associated  with  u"  must  be  significantly  greater  than  the  time 
scale  of  cross-sectional  mixing  due  to  transverse  diffusivity.  A  further 
discussion  of  the  assumptions  made  by  Taylor  and  the  application  of  dis- 
persion coefficients  obtained  in  this  manner  will  be  presented  in  Chapter 
II. 


Having  looked  at  laminar  flows  Taylor  (2)  extended  his  analysis  to 
axisymmetric  turbulent  shear  flow  in  a  pipe.  For  u"  Taylor  used  an  exper- 
imentally verified  velocity  distribution  of  the  form 


u    -u(r) 

-^^ — -  =  Fi(r)  (I-IO) 


where  u    is  the  centerline  velocity,  u+  the  shear  or  friction  velocity, 
max 

and  Fi{r)  an  experimentally  determined  function.  Then  by  assuming  Reynolds 
analogy  which  states  that  the  transfer  of  mass,  heat,  and  momentum  by 
turbulence  are  exactly  analogous,  and  the  same  conditions  required  to 
obtain  (1-7)  he  obtained  a  solution  for  c"  and  numerically  performed  the 
Integration  of  u"c"  shown  in  (1-6).  The  result  was  the  following  expression 
for  the  longitudinal  dispersion  coefficient  in  turbulent  pipe  flow 

E^^  =  10.06  a  u^  (I-ll) 

where  a  is  the  radius  of  the  pipe.  By  assuming  isotropy  of  the  turbulent 
fluctuations  he  was  then  able  to  estimate  the  additional  longitudinal 
transport  due  to  eddy  diffusion  along  the  flow  axis.  Adding  this  to 
(I-ll),  the  corrected  dispersion  coefficient  becomes 

E^  =  10.1  a  u^  (1-12) 

Shortly  after  Taylor's  discovery  and  formulation  of  dispersive  mass 
transport,  Aris  (3)  presented  a  very  elegant  analysis  which  supported 
Taylor's  theory  for  unidirectional  steady  flows  and  rigorously  defined 
the  conditions  under  which  Taylor's  assumptions  were  valid.  Rather  than 
obtaining  solutions  for  u"  and  c"  to  compute  the  dispersive  mass  trans- 
port, Aris  chose  to  apply  a  moment  analysis  of  the  longitudinal  distribution 


of  the  cross-sectional  mean  concentration  of  substance  about  a  point 
moving  with  the  mean  velocity.  The  significance  of  Aris'  analysis  lies 
in  its  generality,  and  the  insight  that  it  yields  for  looking  at  the 
precise  conditions  under  which  an  ordered  correlation  of  u"c"  over  the 
cross-sectional  flow  area  may  be  represented  by  the  implied  stationary 
random  process  of  a  Fickian  flux.  Aris'  development  begins  with  the 
basic  Eulerian  form  of  the  transport  diffusion  equation  for  an  arbitrary 
flow  geometry  and  an  arbitrary  initial  distribution  of  substance.  He 
then  proceeds  to  develop  expressions  for  the  moments  up  through  the  third 
of  the  cross-sectional  mean  concentration  about  5  with  no  assumptions 
regarding  steady  state,  or  the  functional  form  of  the  concentration  c. 
This  generality  allows  him  to  obtain  time  dependent  solutions  which  show 
that  after  an  initial  transient  period,  steady  state  conditions  will  be 
asymptotically  approached  in  which  the  cross-sectional  mean  concentration, 
c,  disperses  about  a  point  traveling  at  the  mean  speed  of  the  flow  in  a 
Gaussian  manner  and  that  the  rate  of  growth  of  the  variance  is  constant 
with  time.  In  other  words  the  critical  assumption  in  Taylor's  analysis 
is  the  one  concerning  steady  state  conditions  with  respect  to  the  moving 
coordinate  system,  i.e.  i^)^  =   0,  which  by  Aris'  analysis  must  be  satis- 
fied if  the  definition  of  the  convective  mass  transport  of  a  substance  as 
a  Fickian  flux  about  the  mean  speed  of  flow  is  to  hold.  The  dimensionless 
time  parameter  governing  the  dispersive  process  is  defined  by  Aris  as 

D^t 
''^iF-  (1-13) 

where  t  represents  time;  D.  is  the  molecular  diffusion  coefficient  in  the 
transverse  direction;  and  a  is  a  characteristic  length  of  the  flow  cross 
section.  The  ratio  a^/D. ,  as  noted  previously,  represents  the  time  scale 


required  for  the  substance  to  mix  over  the  flow  cross  section  due  to  trans- 
verse diffusion.  Thus,  steady  state  conditions  are  approached  for  large 
t,  or  when  the  convective  flow  time  is  large  compared  to  the  time  required 
for  cross-sectional  mixing.  This  supports  Taylor's  comments  regarding  the 
application  of  dispersion  coefficients  to  predict  concentration  distributions. 
In  Chapter  II  it  will  be  shown  how  the  assumption  of  steady  state  condi- 
tions supports  Taylor's  other  assumption  regarding  the  functional  dependence 
of  c  for  both  steady  unidirectional  and  oscillatory  flows.  Aris'  method 
of  analysis  has  been  used  by  other  investigators  to  obtain  expressions 
for  the  dispersion  coefficient  and  for  brevity  later  in  the  text  will  be 
referred  to  as  the  "method  of  moments."  Although  it  remains  a  very   useful 
method  for  determining  expressions  for  longitudinal  dispersion  coefficients 
it  does  not  provide  the  physical  insight  into  the  kinematic  flow  field  and 
cross-sectional  distribution  of  substance  that  Taylor's  method  does. 

Further  support  of  the  general  applicability  of  Taylor's  dispersion 
analysis  to  different  flow  geometries  was  given  by  Elder  (4)  who  applied 
Taylor's  analytical  technique  to  obtain  a  predictive  expression  for  the 
dispersion  coefficient  in  two  dimensional  unidirectional  shear  flows  in 
infinitely  wide  open  channels.  In  his  analysis  Elder  obtains  the  following 
expression  for  the  longitudinal  dispersion  coefficient: 

E^  -  ^  \    u"  \   i   \    u"  dz'dz'dz'  (1-13) 

where  z'  is  a  variable  of  integration;  z  is  the  vertical  coordinate;  h  is  the 
depth  of  flow;  and  e  is  either  the  vertical  molecular  diffusivity  or  the 
vertical  eddy  diffusivity  depending  upon  whether  the  flow  is  laminar  or 
turbulent.  This  expression  indicates,  in  a  manner  similar  to  (1-8), 


that  in  either  laminar  or  turbulent  unidirectional  shear  flow  the  longi- 
tudinal dispersion  varies  directly  with  the  time  scale  of  cross-sectional 
mixing.  To  obtain  an  expression  for  u"  Elder  assumed  the  velocity  distri- 
bution to  be  logarithmic  which  upon  substitution  into  (1-13)  and  correcting 
for  longitudinal  eddy  diffusion  as  done  by  Taylor,  yields 

El  =  5.93  u^h  (1-14) 

where  u*  is  the  shear  velocity. 

The  importance  of  the  flow  cross  section  geometry  on  the  velocity 
shear  profile,  transverse  diffusivity,  and  ultimately  the  longitudinal 
dispersion  is  demonstrated  by  (1-12)  and  (1-14)  which  if  expressed 
in  terms  of  the  hydraulic  radius,  Rn,  for  the  respective  flows  yield  values 
of  E. /R„u^  of  20.2  for  the  pipe  and  5.9  for  the  infinitely  wide  channel. 
Bowden  (5)  developed  expressions  for  E,  using  several  different  velocity 
and  eddy  diffusivity  distributions  with  depth  and  showed  that  E./u^h 
ranged  from  5.9  to  25  for  the  cases  analyzed.  Thus,  the  particular  velocity 
shear  profile  used  along  with  the  cross-sectional  eddy  diffusivity  have 
a  significant  effect  on  the  magnitude  of  the  predicted  longitudinal 
dispersion  coefficient. 

While  the  works  of  Taylor,  Aris,  and  Elder  proved  the  existence  of 
an  effective  dispersive  mass  transport  mechanism  and  there  was  reasonable 
agreement  between  laboratory  experimental  data  for  pipe  and  infinitely 
wide  open  channel  flow  geometries,  other  investigators  were  finding  that 
observed  values  of  the  longitudinal  dispersion  coefficient  for  natural 
rivers  and  streams  were  considerably  higher  than  those  predicted  by  theory. 
As  reported  by  Fischer  (6)  observed  values  for  E.  in  natural  watercourses 
ranged  from  50  to  700  hu^  as  compared  with  the  5.9  to  25  hu^  range  predicted 
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by  Elder  and  Bowden.  To  account  for  this  increased  effect  Fischer  concluded 
that  the  dominant  mechanism  in  longitudinal  dispersion  was  the  interaction 
of  the  velocity  shear  profile  and  turbulent  mixing  time  scales  across  the 
finite  width  of  the  channel  as  opposed  to  the  vertical  variations  treated 
by  the  previous  investigators  for  infinitely  wide  channels.  The  effect 
of  asymmetrical  flow  cross  sections  on  longitudinal  dispersion  was  first 
examined  by  Aris  (3)  who  showed  that  the  dispersion  in  a  circular  tube  is 
less  than  that  in  an  elliptical  one  of  the  same  area.  This  in  addition  to 
the  fact  that  all  real  watercourses  have  a  finite  width  and  that  width- 
to-depth  ratios  for  natural  channels  are  usually  significantly  greater 
than  unity  makes  Fischer's  hypothesis  very  reasonable.  For  wide  channels 

2 

the  transverse  mixing  time,  i  /K  ,   over  some  characteristic  i    would  tend 
to  be  greater  than  the  corresponding  vertical  mixing  time  scale,  h^/K  . 
Fischer  thus  argued  that  the  velocity  shear  profile  across  the  channel 
would  have  a  greater  longitudinal  dispersive  effect  since  an  increased 


mixing  time  would  produce  larger  values  of  the  spatial  correlation  u"c". 
Using  this  as  the  basis  for  his  analysis  Fischer  applied  Taylor's  technique 
to  a  rectangular  coordinate  system  and  assumed  that  vertical  variations 
in  c"  were  negligible  compared  to  the  transverse  variations  to  obtain 

""!  =  ¥  '^f">  <'-'^' 

where  y  is  the  coordinate  axis  across  the  channel. 

Fischer  then  integrated  (1-15)  over  the  cross-sectional  area  of  the 
channel  to  obtain  an  expression  for  c"  which  he  then  correlated  with  u" 
for  the  dispersive  mass  transport  and  using  Taylor's  definition  of  a 
Fickian  flux  obtained  a  longitudinal  dispersion  coefficient  of 
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^L  =  -l  \  ^"^^^''^  \   rhiyy^^  \   ''"^y^'^y         (^-''^^ 


where  the  depth,  h,  is  a  function  of  y;  K  is  the  transverse  eddy  diffusi- 
vity;  w  is  the  width  of  the  channel;  and  q"(y)  is  the  depth  integrated 
flow  defined  by 

h(y) 

q"(y)  =  ^    u"(y,z)  dz  (1-17) 

0 

To  obtain  a  more  useful  expression  for  E.  Fischer  defines  a  Lagrangian 
time  scale  for  cross-sectional  mixing  which  he  obtains  through  an  extension 
of  Taylor's  work  on  diffusion  by  continuous  movements  (7),  and  relates  it 
to  the  Eulerian  time  scale  for  cross-sectional  mixing  previously  discussed. 
Then,  using  Elder's  (4)  experimental  determination  that  the  transverse 
eddy  diffusivity  in  an  infinitely  wide  channel  with  a  steady  unidirectional 
flow  can  be  expressed  as 

K  =  0.23  hu*  (1-18) 

Fischer  arrives  at  an  alternate  expression  for  E.  given  by 

where  Z     in  a  natural  channel  is  the  distance  from  the  point  of  maximum 
s 

surface  velocity  to  the  most  distant  bank.  To  verify  his  theory  Fischer 
conducted  a  detailed  set  of  laboratory  experiments.  Using  width-to-depth 
ratios  ranging  from  9  to  15.7  Fischer  obtained  good  agreement  between 
predicted  and  observed  values  of  E.  using  (1-16).  However,  (1-19)  over- 
predicted  in  each  case,  with  a  maximum  error  of  75  per  cent. 

Fischer  (8)  also  applied  his  analysis  to  data  obtained  by  other  investigators 
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in  6  natural  watercourses  including  Copper  Creek,  Va.  and  Clinch  River, 
Tenn.  (2  locations  each),  Powell  River,  Tenn.,  and  Coachella  Canal,  Cali- 
fornia. Width-to-depth  ratios  for  these  sites  range  from  15  to  62  with 
most  being  less  than  40.  To  carry  out  his  analysis  Fischer  numerically 
integrated  (1-16)  using  measured  data  and  found  that  in  uniform  channels 
the  predicted  values  for  E^  were  within  30  per  cent  of  the  observed 
values  whereas  in  non-uniform  channels  the  predicted  value  varied  from 
the  observed  by  as  much  as  a  factor  of  4. 

Thus,  while  Fischer's  hypothesis  regarding  the  dominance  of  lateral 
effects  on  longitudinal  dispersion  in  unidirectional  flows  yields  reason- 
able results  his  analysis  provides  very  little  insight  into  the  physics 
of  the  three  dimensional  problem.  Both  (1-16)  and  (1-19)  make  no  attempt 
to  describe  the  velocity  distribution  over  the  flow  cross  section  thereby 
requiring  detailed  measurements  for  their  application.  Moreover,  by  neg- 
lecting the  vertical  shear  effects  it  remains  impossible  to  obtain  from 
his  analysis  the  conditions  for  which  either  the  vertical  or  lateral 
effects  would  be  the  dominant  mechanism  for  longitudinal  dispersion. 
I-B.2  Oscillatory  Flow 

The  analysis  of  dispersive  mass  transport  in  oscillatory  flow  is 
more  complex  than  the  corresponding  problem  in  steady  unidirectional  flow 
because  of  the  unsteady  nature  of  the  governing  equations  of  motion  and 
transport  diffusion.  Since,  however,  many  important  water  quality  problems 
involve  the  predicted  distribution  and  transport  of  substances  within 
tidal  waters  the  analysis  of  dispersive  processes  in  oscillatory  flows 
has  attracted  considerable  interest.  There  is  in  fact  a  reasonable  doubt 
that  one  is  justified  in  representing  the  convective  mass  transport  u"c" 
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as  a  Fickian  flux  in  an  oscillatory  flow.  However,  it  will  be  shown  in 
Chapter  II  that  under  the  assumptions  of  periodicity  and  uniform  flow  the 
analyses  of  Taylor  and  Aris  apply. 

Bowden  (5)  was  the  first  investigator  to  look  at  longitudinal  dis- 
persive mass  transport  in  oscillatory  flow.  Bowden  generally  followed 
Taylor's  original  analytical  technique  and  began  by  assuming  a  monochro- 
matic flow  in  an  infinitely  wide  channel  in  which  u"  can  be  expressed  as 

u"'(z,t)  =  Ai(z)  cos  at  +  Bi(z)  sin  at  (1-20) 
He  likewise  assumed  that  c"  is  simple  harmonic  in  time  and  of  the  form 

c"(z,t)  =  Pi(z)  cos  at  +  Qi(z)  sin  at  (1-21) 
conditional  upon  satisfying  the  transport  diffusion  equation  expressed  as 

This  is  a  very  proper  formulation  of  the  problem,  allowing  for  the 
periodic  variation  of  u"  and  c"  with  time  and  the  existence  of  temporal 
phase  shifts  for  u"  and  c"  as  a  function  of  position  within  the  water 
column.  Unfortunately,  Bowden  proceeded  by  means  of  his  assumptions  to 

reduce  the  problem  to  one  which  applies  only  to  the  case  for  which  the 

3c" 
period  of  oscillation  is  infinitely  long,  {j^)r  =   0,  and  for  which  the 

temporal  phase  shift  as  a  function  of  position  within  the  water  column 
has  been  eliminated.  For  these  conditions  Bowden  found  that  the  longitu- 
dinal dispersion  coefficient  in  an  oscillatory  flow  is  one-half  the  value 
of  the  same  coefficient  for  a  corresponding  unidirectional  flow  having  the 
same  surface  velocity,  the  same  shear  profile,  and  the  same  vertical  eddy 
diffusivity. 
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Holley  and  Harleman  (9)  treated  the  case  of  oscillatory  turbulent  pipe 
flow  and  obtained  an  expression  for  the  longitudinal  oscillatory  flow  dis- 
persion coefficient,  E  ,  using  Aris'  method  of  moments.  In  their  analysis 
they  assumed  that  at  each  instant  in  time  the  lateral  distribution  of 
velocity  and  eddy  diffusivity  are  the  same  as  they  would  be  if  the  flow 
were  steady  with  the  same  cross-sectional  mean  velocity.  This  assumption 
effectively  removes  the  oscillatory  nature  of  the  flow  and  reduces  the 
problem  to  essentially  the  same  one  solved  by  Aris  (3)  and  Taylor  (2). 
Thus,  it  is  not  surprising  that  Holley  and  Harleman  predicted  a  longitu- 
dinal dispersion  coefficient  of 

E  (t)  =  10.1a  u*  (t)  (1-23) 

where,  as  before,  a  is  the  radius  of  the  pipe,  and  u^  (t)  is  the  friction 
velocity  now  periodic  in  time.  Predicted  values  of  E^  were  compared 
with  observations  made  from  a  series  of  laboratory  experiments  in  which 
an  oscillatory  flow  was  produced  in  a  1-1/2  inch  pipe  by  a  piston  genera- 
tor. Superimposed  upon  the  oscillating  flow  was  a  small  unidirectional 
flow  to  simulate  the  net  fresh  water  discharge  of  an  estuarine  type  flow. 
The  experimental  data,  however,  showed  that  due  to  the  small  magnitude 
of  the  unidirectional  flow  compared  to  the  velocity  amplitude  of  the  oscilla- 
tory flow,  its  effect  on  longitudinal  dispersion  could  be  neglected.  The 
data  also  showed  that  (1-23)  generally  underpredicted  the  dispersion 
coefficient  with  better  agreement  between  predicted  and  observed  values 
coming  with  higher  Reynolds  number  flows.  The  authors  attributed  this 
behavior  to  laminar  sublayer  effects. 

Okubo  (10)  obtained  a  less  restrictive  solution  which  for  the  first 
time  demonstrated  that  in  a  two  dimensional  oscillatory  shear  flow  the 


15 


period  of  oscillation  and  the  time  scale  of  vertical  mixing  played  an 
important  role  in  the  dispersive  process.  Okubo  confined  his  analysis 
to  vertically  bounded  and  unbounded  flows  with  no  transverse  shear.  Only 
the  former  will  be  discussed.  He  assumed  a  linear,  monochromatic  flow 
superimposed  upon  a  linear  unidirectional  flow,  or 

u(z,t)  =  Uo  (1-|)  +  Vo  (1-|)  sin  at  (1-24) 

where  Uo  and  Vo  are  the  surface  amplitudes  of  the  unidirectional  and 
oscillatory  flow,  and  a   is  the  angular  frequency  of  oscillation.  Okubo 
then  assumed  a  constant  eddy  diffusivity  both  in  the  direction  of  flow 
and  over  the  depth,  and  zero  diffusive  flux  conditions  at  the  surface  and 
bottom.  Using  the  method  of  moments  as  applied  by  Aris  and  the  assumptions 
discussed  he  obtained  an  expression  for  the  longitudinal  dispersion  coeffi- 


cient, E  ,  functionally  expressed  as 


h^ 


E,  =  X  (Uo,  Vo,  h,  T,  ^)  (1-25) 

X  ^2 

1,2 

where  T  is  the  period  of  oscillation,  and  ^  the  time  scale  of  vertical 
mixing.  To  interpret  his  results  Okubo  considered  two  extreme  cases. 

(1)  S-  »  T  -  Neglecting  longitudinal  eddy  diffusion  his  solution  becomes 
■^z 

ui    .2  n    1% 

unidirectional         oscillatory 
From  this  Okubo  correctly  concluded  that  for  situations  in  which  the  verti- 
cal mixing  time  is  large  compared  with  the  period  of  oscillation  the  con- 
tribution to  the  longitudinal  dispersive  mass  transport  due  to  the  oscilla- 
tory flow  is  small  compared  with  the  contribution  due  to  the  unidirectional 
flow.  However,  (1-26)  also  suggests  that  the  functional  behavior  of  dispersion 
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in  an  oscillatory  flow  is  significantly  different  from  the  same  process  in 
unidirectional  flow.  One  of  the  principal  findings  of  the  present  work 
concerns  this  point  and  will  be  treated  in  depth  in  Chapters  III  and  IV. 


(2)  T  »  p-  -  Again  neglecting  longitudinal  eddy  diffusion  Okubo's  result 


may  be  written  as 

E  =  !!!_  (hi)  +  !L  (hi)  (1-27) 

S   120  ^K^^   236  ^K^^  ^    '' 

This  case  is  comparable  to  Bowden's  analysis  for  a  wery   long  period  of 
oscillation  and  supports  his  finding  that  E^  ->  1/2  E|_  for  large  T. 

Holley  et  al .  (11)  successfully  applied  Taylor's  original  method  of 
analysis  to  obtain  an  analytical  expression  for  <E^>j  in  a  two  dimensional 
periodic,  uniform  shear  flow.  For  their  analysis  they  assumed  the  spatial 
variation  component  of  velocity  to  be 

u"  (z.t)  =  a  z  sin  at  (1-28) 

where  ct  is  a  constant.  This  is  identical  to  the  linear  oscillatory  flow 
profile  assumed  by  Okubo.  This  expression  for  u"  is  then  used  to  force 
the  following  form  of  the  transport  diffusion  equation  to  obtain  a  solution 
for  c"  (z.t) 

/8c'\    ,,  3^c"  _  „„  8c  /T  29) 

The  solutions  obtained  for  u"  and  c"  are  then  used  to  determine  the  depth 
integrated  time  mean  convective  mass  transport  with  respect  to  the  depth 
mean  velocity.  The  resulting  expression  for  the  mean  longitudinal  dispersion 
coefficient  over  the  period  of  oscillation  is 

<E>  =5CTiJil  ^   _! (i_3o) 

^  ^   ^  S    n=l   (2n-l)M[|  (2n-l)^  J'r+^} 


17 


where  T'  =  -r4.  Hoi  ley  et  al .  then  proceed  to  relate  <E^>j  in  an  oscillating 

flow  to  the  dispersive  process  in  an  oscillating  flow  of  infinite  period 

with  the  same  velocity  profile  as  given  by  (1-28).  Solving  for  c"  using 

the  steady  state  form  of  (1-29)  and  proceeding  in  the  usual  manner  they 

obtain  an  expression  for  <E  >t  which  they  call  E^.  It  is  shown  that  E^ 

is  one-half  the  corresponding  value  for  E,  in  agreement  with  Bowden  (5), 

and  is  subsequently  used  to  compare  <E^>-p  with  the  corresponding  steady 

unidirectional  flow  coefficient.  Apparently  the  investigators'  primary  goal 

was  not  to  describe  longitudinal  dispersion  in  an  oscillatory  flow  as  a 

separate  distinct  physical  process  but  rather  to  describe  it  in  terms  of  a 

"corresponding"  process  in  a  unidirectional  flow.  To  do  this  the  solution 

for  <E  >j  given  by  (1-30)  was  normalized  by  E^  and  the  result  plotted  against 

TK 
the  non-dimensional  time,  T'  =  -p^,  suggested  by  Okubo's  (10)  results. 

Numerical  simulations  were  also  carried  out  for  other  velocity  profiles 

in  which  the  velocity  was  again  assumed  to  be  temporally  in  phase  over 

depth  during  the  period  of  oscillation.  Based  upon  their  analytical  and 

numerical  results  Holley  et  al .  concluded  that: 

(1)  T'>1  -  The  longitudinal  dispersion  coefficient  in  an  oscillatory 
flow  is  independent  of  T'  and  equal  to  one-half  the  value  of  the  dispersion 
coefficient  in  a  "corresponding"  unidirectional  flow. 

(2)  T'<0.1  -  The  longitudinal  dispersion  coefficient  in  an  oscilla- 
tory flow  rapidly  becomes  insignificant  compared  to  the  dispersion  coefficient 
in  a  "corresponding"  unidirectional  flow  and  is  functionally  described  by 

<E  >^  .  T'2  (1-31) 

X  T 

These  results  have  been  subsequently  applied  by  Fischer  (12,  13),  Awaya 

(14),  and  Fischer  and  Holley  (15)  in  analyses  of  dispersion  in  oscillatory 
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flows.  Fischer  (12)  attributes  the  functional  behavior  of  oscillatory  flow 
dispersion  to  a  temporal  phase  shift  between  u"  and  c"  as  a  function  of  T'. 
However,  it  will  be  shown  in  the  present  work  that  while  these  conclusions 
are  applicable  for  the  non-dimensional  parameter  — f — ,  they  are  incorrect 
and  misleading  in  their  description  of  the  longitudinal  dispersive  mass 
transport  process  in  oscillatory  flow. 

Fukuoka  (16)  obtained  analytical  expressions  for  E^  from  several  cases 
of  two  dimensional  and  axisymmetric  oscillatory  shear  flows  using  the  method 
of  moments.  For  his  analysis  he  assumed  a  velocity  of  the  general  form 

u  =  u  (z)  +  u^  (z)  sin  at  (1-32) 

similar  to  Okubo  (10),  implying  once  again  no  temporal  phase  shift  of  the 
velocity  over  the  flow  cross  section.  In  specifying  the  forms  of  u^  (z) 
and  Uo.  (z)^  Fukuoka  assumed  that  the  spatial  variation  at  each  instant  of 
time  is  the  same  as  that  for  an  equivalent  unidirectional  flow.  Specific 
cases  treated  in  depth  are: 

(1)  linear  velocity  profile  over  depth  with  constant  K^ 

(2)  parabolic  velocity  profile  for  axisymmetric  flow  with  constant 


h 


sin  at 


(3)  linear  velocity  profile  over  depth  with  K^  =  Kc 

In  interpreting  his  results  Fukuoka  follows  Holley  et  al .  by  plotting 

•=^  vs.  T'  and  arrives  at  the  same  conclusion  regarding  the  functional  be- 
havior  of  dispersion  in  oscillatory  flows.  For  Case  3  he  shows  analytically 
that  the  dispersion  coefficient  corresponding  to  a  time  varying  eddy  diffu- 
sivity  differs  by  a  factor  of  S/tt^  from  one  that  is  considered  to  be 
constant,  provided  the  constant  K^  is  |  Ko.  This  small  effect  of  a  periodic 
eddy  diffusivity  supports  results  obtained  numerically  by  Holley  et  al . 
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It  is  interesting  to  note  that  none  of  the  analyses  dealing  with 
oscillatory  flow  discussed  thus  far  have  attempted  to  solve  analytically 
the  velocity  distribution  over  the  flow  cross  section.  Moreover,  it  has 
been  further  assumed  that  the  velocity  profile  at  each  instant  in  time  is 
the  same  as  an  equivalent  unidirectional  flow  profile  and  that  no  temporal 
phase  shift  as  a  function  of  spatial  position  exists.  These  assumptions 
are  particularly  unrealistic  when  one  considers  the  definition  of  dispersive 
mass  transport  as  defined  by  Taylor,  i.e-JL  u"c"  dA.  By  restricting  the 
form  of  u"  in  this  manner  one  is  also  restricting  the  form  of  c"  since  by 
(1-29)  it  is  seen  that  the  distribution  of  substance  is  in  part  forced  by 
the  convective  flux.  Thus,  any  artificiality  regarding  the  assumed  form 
of  u"  appears  in  the  dispersion  coefficient  through  both  u"  and  c".  The 
significantly  different  nature  of  a  truly  oscillatory  flow  as  compared  with 
a  unidirectional  one  is  illustrated  by Schlichting  (17)  in  his  treatment 
of  the  fluid  motion  above  an  oscillating  flat  plate,  and  Lamb  (18)  in  his 
treatment  on  the  effect  of  bottom  friction  on  long  waves  in  two  dimensional 
shear  flow.  Although  the  forcing  functions  in  these  two  cases  are  different 
(oscillating  boundary  shear  in  Schlichting  as  opposed  to  an  oscillating 
pressure  gradient  in  Lamb),  the  result  of  a  shear  wave  propagating  through 
the  fluid  in  a  direction  normal  to  the  axis  of  flow  is  similar,  causing  a 
continuously  changing  phase  shift  of  the  flow  over  the  cross-sectional  area. 

Awaya  (14),  Segall  (19),  and  Segall  and  Gidlund  (20)  have  obtained 
analytical  expressions  for  the  longitudinal  dispersion  coefficient  in 
oscillatory  flow  situations  using  realistic  velocities  obtained  directly 
from  solutions  to  the  appropriate  forms  of  the  equations  of  motion.  Awaya 
addresses  the  case  of  an  oscillatory  laminar  flow  in  a  circular  cross  section 
for  which  he  obtains  an  expression  for  <E  >j  using  Taylor's  technique,  and 
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normalizes  it  by  Ea,  after  Holley  et  al.  His  results  are  similar  to  those 
of  Holley  et  al .  and  his  interpretation  the  same.  Segall  and  Segall  and 
Gidlund  use  a  simplified  version  of  Lamb's  previously  mentioned  solution 
for  u  corresponding  to  small  values  of  the  vertical  eddy  viscosity,  e^, 
which  they  apply  to  the  transport  diffusion  equation  and  solve  using  the 
method  of  moments.  In  both  of  these  analyses  the  interpretation  of  the 
functional  behavior  of  the  dispersive  process  in  oscillatory  flows  is  the 
same  as  that  of  previous  investigators  although  the  matter  is  not  pursued 
in  any  detail . 
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I-C  Description  and  Scope  of 
the  Present  Work 

In  the  chapters  that  follow  four  sets  of  boundary  value  problems  will 
be  solved  yielding  analytical  solutions  for  u",  c",  E|_,  and  <E^>j  for  the 
following  cases: 

(1)  Steady  unidirectional  flow  in  an  infinitely  wide  rectangular 
channel ,  e  and  K  constant 

(2)  Monochromatic  periodic  flow  in  an  infinitely  wide  rectangular 
channel ,  e  and  K  constant 

(3)  Steady  unidirectional  flow  in  a  rectangular  channel  of  width 
w;  e^,  e  ,  K^.  and  K  constant 

(4)  Monochromatic  periodic  flow  in  a  rectangular  channel  of  width 
w;  e^,  e  ,   K^,  and  K  constant. 

These  solutions  deal  exclusively  with  dispersive  mass  transport  in  homo- 
geneous fluids  and  uniform  or  nearly  uniform  flows. 

In  Chapter  II  a  discussion  is  presented  of  the  conditions  under  which 
Taylor's  definition  of  dispersive  mass  transport  is  valid  for  unidirectional 
and  oscillatory  flows.  Chapter  III  treats  Cases  1  and  2  above  yielding 
solutions  for  E.  and  <E  >-r  in  two  dimensional  shear  flows  while  Chapter  IV 
extends  the  analysis  to  three  dimensional  shear  flows  through  Cases  3  and  4. 
In  Chapter  V  predicted  values  of  E^  and  <E^>y  using  expressions  developed 
in  Chapters  III  and  IV  are  compared  with  available  data.  Chapter  VI  presents 
a  summary  of  the  present  work  and  some  conclusions  drawn  from  it.  The  mathe- 
matical details  associated  with  portions  of  the  boundary  value  problem 
solutions  are  reserved  for  Appendices  A,  B,  and  C. 


CHAPTER  II 

CONSIDERATIONS  FOR  APPLICATION  OF  TAYLOR'S 

DISPERSION  ANALYSIS  TO  FREE  SURFACE  FLOWS 

II-A  General 


In  order  to  provide  a  solid  base  upon  which  the  analyses  presented  in 
Chapters  III  and  IV  are  built,  it  is  necessary  to  set  forth  at  this  juncture 
a  detailed  development  of  Taylor's  dispersion  analysis  as  it  applies  to 
both  unidirectional  and  oscillatory  free  surface  flows.  By  doing  this  it 
is  hoped  that  all  assumptions  made  in  obtaining  solutions  for  u",  c",  E 
and  E.  ,  and  the  conditions  under  which  a  dispersion  coefficient  may  be  defined 
and  applied  consistent  with  Aris  (3)  and  Taylor  (1,2),  will  be  clearly 
understood. 

Consider  a  three  dimensional  shear  flow,  either  unidirectional  or 
oscillatory,  in  a  rectangular  channel  of  width  w  and  depth  h.  Let  the 
coordinate  system  be  selected  as  shown  in  Figure  1  with  the  x  coordinate 
taken  along  the  longitudinal  axis  of  the  channel  denoting  the  direction 
of  the 


£. 


J^ 


u  (y.z.t) 


y«-w/2- 


z  =  0 


^ 


Tn:^ 


y  =  w/2 


FIGURE  1  DEFINITION  SKETCH 
velocity  component,  u(y,z,t);  the  y  coordinate  denoting  lateral  position 
across  the  channel;  and  the  z  coordinate  acting  positive  upward  from  the 
free  surface.  To  begin  the  analysis  it  shall  be  assumed  that: 
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(1)  Longitudinal  mass  transport  due  to  eddy  diffusion  along  the  flow 
axis  is  small  compared  to  convective  transport  mechanisms. 

(2)  Eddy  diffusivities  K  and  K  can  be  represented  by  their  mean 
values  over  the  flow  cross  section, 

(3)  Uniform  flow  conditions  along  the  principal  axis  of  the  channel 
exist  as  indicated  by  the  functional  notation  for  u. 

Applying  these  assumptions  to  the  channel  shown  in  Figure  1  the 
general  form  of  the  transport  diffusion  equation  for  a  conservative  sub- 
stance may  be  written  as 

|£+  u|f  =  "k  ^+  "k  ^  (ll-l) 

3t    3x   y  3y2   z  g^2 

After  Taylor  (1,2),  transform  (II-l)  to  a  coordinate  system  traveling 

with  the  cross-sectional  mean  velocity  using 

u  =  u  +  u"  (II-2) 

C  =  x  -  ut  (II-3) 

so  that  the  functional  transformation  of  the  concentration  c  is  given  by 
c(x,y,z,t)  ^  c(C(x,t),y,z,t). 

The  transformation  relationships  are  then 


3c  _  / 3Cv   .  9c  8i 
8t  "  ht'^       dE,  dt 


3x   3C   3x 


—  and  ^Y  unchanged 
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Using  (II-3), 


9c  _  ,3cx    -  ac  (II-4) 


3£  _  9£ 
3x  "  3^ 


(II-5) 


Equations  (II-2),  (II-4),  and  (II-5)  are  then  substituted  into  (II-l)  to 
yield 

(^)_  +  u"^='k^+~K^  (II-6) 

Next,  decompose  c  into  its  cross-sectional  mean  and  variation  components 
by  the  expression 

c  =  c  +  c"  (11-7) 

and  substitute  into  (II-6).  For  uniform  flow  this  produces 

/3C\    ^3c"\    Z     3  c"  z    3  c"  _   3  (u"c")    II  3c      /yr  n^ 

which  is  the  form  of  the  transport  diffusion  equation  as  seen  from  a 
coordinate  system  traveling  with  the  mean  cross-sectional  velocity. 

Proceeding  with  the  analysis,  (II-8)  is  integrated  over  the  cross- 
sectional  area, 

Ji  w/2    _        h  w/2  h  w/2   ^ 

0  -w/2  0  -w/2  0  -w/2 

h  w/2  h  w/2  h  w/2 

0  -w/2  0  -w/2  0  -w/2 

It  is  next  assumed  that: 

(1)  The  time  rate  of  change  of  the  depth  is  small. 
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(2)  There  is  no  turbulent  diffusion  of  substance  across  the  enclosing 

boundaries  of  the  flow  cross  section. 

To  apply  the  condition  of  zero  diffusive  flux  across  the  boundaries 
of  the  channel  it  is  helpful  to  use  Green's  Theorem  in  the  plane  which  for 
the  third  and  fourth  terms  on  the  left  hand  side  of  (II-9)  states 

A  u 

where  u  is  the  perimeter  of  the  area  A.  Since  by  assumption  (2)  it  is 
required  that  there  be  no  diffusive  flux  across  u,  the  left  hand  side  of 
(11-10)  is  equal  to  zero.  Using  this  result,  assumption  (1),  and  noting 
that 


u"  dA  =  \\c"  dA  =  0 
A       A 
by  definition,  (II-9)  is  reduced  to 


As  noted  in  I.B.I  the  form  of  (11-11)  induced  Taylor  to  define  the  spa- 


tially averaged  convective  flux  u"c"  as  a  Fickian  flux  thus  transforming 
(11-11)  to  the  one  dimensional  heat  equation  form 

,8Cx  ^  ^  3fc  (11-12) 

where. 


•  >i_  II 


^L  =  711,  '"-^" 
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Referring  again  to  I.B.I,  Aris  showed  that  (11-12)  and  (11-13)  are 
valid  formulations  if  the  following  necessary  and  sufficient  conditions 
are  satisfied:  (1)  steady  uniform  flow;  and  (2)  steady  state  concentration 
with  respect  to  a  point  traveling  at  the  mean  velocity,  i.e.  (rv-),.  =  0. 
In  the  sections  that  follow  it  will  be  shown  that  these  conditions  are 
necessary  and  sufficient  not  only  for  unidirectional  flows  but  for  oscilla- 
tory flows  as  well  provided  that  steady  state  is  defined  in  the  periodic 
sense.  It  will  also  be  shown  that  Taylor's  assumptions  regarding  the 
functional  dependence  of  c  and  c"  on  E,,   y,  and  z  as  stated  in  I.B.I 
follow  directly  from  the  more  general  requirements  of  steady  state  u  and 
c,  and  uniform  flow. 
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II-B  Unidirectional  Flow 

To  obtain  an  expression  for  E.  as  defined  by  (11-13)  it  is  necessary 
first  to  have  solutions  for  u"  and  c".  Ideally  a  solution  for  u"(y,z)  is 
obtained  from  the  governing  equation  of  motion  in  the  x  direction.  This 
solution  is  then  used  to  force  a  simplified  form  of  (II-8)  to  obtain  a 
solution  for  c"  which  is  then  correlated  with  u"  according  to  (11-13)  for 
the  determination  of  the  longitudinal  dispersion  coefficient.  As  stated 
in  I.B.I  Taylor  reduced  an  expression  similar  to  (II-8)  to  a  solvable 
form  by  assuming: 

(1)  c"  is  not  a  function  of  E,. 

3c 

(2)  -r^  is  constant. 

OS 

(3)  Steady  state  conditions  with  respect  to  the  moving  coordinate 
system  exist,  i.e.  (||-)^  =  0. 

Aris,  however,  showed  that  assumption  (3)  along  with  the  requirements 
of  steady  uniform  flow  were  necessary  and  sufficient  for  the  definition 
of  a  dispersion  coefficient  as  stated  by  (11-13),  and  it  will  now  be 
demonstrated  that  assumptions  (1)  and  (2)  follow  from  these  conditions. 

For  steady  state  conditions  to  exist  (11-11)  reduces  to 

A 

Assuming  the  channel  to  be  prismatic,  the  order  of  differentiation  and 
integration  may  be  reversed  yielding 

0  =  U   1^  (u"c")  dA 
which  for  uniform  flow  becomes 
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0  =   \\  u"  ll^dA  (11-14) 


To  examine  the  functional  dependence  of  c",  let 
u"  =  92  {y*2) 

If"  =  f2  (y.z.U 
Equation  (11-14)  may  then  be  written  as 

0  =  U  92  (y,z)-f2  (y,z,0  dA  (11-15) 

which  upon  integration  yields 

G(U  =  0  (11-16) 

From  (11-16)  it  may  be  inferred  that  ||-  =  0  or  c"  is  not  a  function  of 
5  thereby  justifying  Taylor's  first  assumption  restated  as 

c(x.y,z)  =  E(x)  +  c"(y,z)  (11-17) 

The  above  result  is  then  used  to  substantiate  Taylor's  second  assump- 
tion that  ||-  is  constant.  To  do  this  (II-8)  is  differentiated  with 
respect  to  E,.     For  steady  state  conditions  and  uniform  flow  the  result 
is 

u"|^=0  (11-18) 

OS. 

which  upon  integration  yields 


|^=  constant  (H-IS) 

The  application  of  (11-17)  and  (11-19)  reduces  (II-8)  to 
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where  ^  is  constant. 

Thus  it  is  shown  that  the  definition  of  the  longitudinal  dispersion 
coefficient  as  stated  by  (11-13),  which  implies  the  use  of  solutions  for 
u"  and  c"  obtained  from  the  governing  equation  of  motion  and  (11-20) 
respectively,  is  valid  for  free  surface  flows  provided  the  necessary  and 
sufficient  conditions  of  uniform  steady  flow  and  steady  state  concentration 
distributions  with  respect  to  a  coordinate  system  traveling  with  the  mean 
velocity  are  satisfied.  Taylor's  assumptions  as  stated  by  (11-17)  and 
(11-19)  follow  from  these  conditions. 
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II-C  Oscillatory  Flow 

The  same  analysis  carried  out  in  II-B  for  steady  unidirectional  flow 
can  be  extended  to  the  oscillatory  flow  case  if  it  is  assumed  that  steady 
state  conditions  in  the  oscillatory  sense  can  be  satisfied  by  requiring 
both  the  velocity,  u,  and  the  concentration,  c,  to  be  temporally  periodic. 
This,  of  course,  implies  that  u,  c,  u",  and  c"  are  also  temporally  periodic. 

To  show  that  Taylor's  same  assumptions  apply  to  dispersion  in  oscilla- 
tory flows  it  is  again  necessary  only  to  require  that  the  flow  be  uniform 
and  that  the  velocity  and  concentration  be  steady  state  functions  in  the 
oscillatory  sense,  i .e.,  periodic.  If  this  can  be  shown  to  be  true  then 
it  is  reasonable  to  assume  that  Aris'  argument  could  be  extended  to  periodic 
functions  thereby  justifying  the  definition  of  a  longitudinal  dispersion 
coefficient  in  oscillatory  flow  using  Taylor's  method. 

The  analysis  is  begun  by  averaging  (11-11)  over  one  full  period  of 
oscillation  which  for  a  periodic  c  yields 

<-  ^  (UVO>,  =  0  (11-21) 

where  C  is  now  defined  by 

^  =  X  -   \  D  (f)  df  (11-22) 

-^0 

Next,  it  is  assumed  once  again  that  the  time  rate  of  change  in  depth  is 

small  so  that  the  implied  order  of  integration  and  differentiation  in 

(11-21)  may  be  reversed 

ph       ^w/2       pT 
^        \  \  \        u"c"  dt  dy  dz  =  0  (11-23) 

-' 0       J-w/2     Jq 

Since  it  is  required  that  the  flow  be  uniform  and  that  u  and  c  be 
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temporally  periodic,  assume 

u"  =  gi(y.z)  cos  at 

c"  =  fi(y,z,c)  cos(at  -ai) 
where  ai  is  a  constant.  It  then  follows  that  the  temporal  integration 
of  u"c"  is  of  the  form 

\  u"c"dt  =  B(y,z,5)  cos  ai  (11-24) 

Substituting  (11-24)  into  (11-23) 
ph   ^w/2 

3l 


ph   pW/2 
|-  \    \    B(y,z,c)  cos  ai  dy  dz  =  0 


0    -w/2 
which  upon  integration  reduces  to  the  form 

cos  a,   •  f  (^^  =  0  (11-25) 

From  this  it  is  seen  that  X(^)  is  a  constant  and  thus  c"  is  not  a 
function  of  C-  This  is  the  same  result  obtained  in  II-B  and  is  consistent 
with  Taylor's  first  assumption.  However,  for  the  oscillatory  case  there 
is  an  additional  distinction  to  be  made.  Recognizing  that  uniform  flow 
is  required  and  c"  is  not  a  function  of  C,  (11-11)  reduces  to 

and  (1 1-8)  becomes 

,dc'\         Z     9V:   K  ^^=  -  u"  ^  (11-27) 

^^C  ■  >  3y2   "^z  3Z2    "^  9C  ^^^  ^  ' 

Differentiating  (11-27)  with  respect  to  C  yields 

u"  lie  =  0 
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or,  ||-  =  constant,  which  is  again  consistent  with  Taylor's  assumption  for 

3c 

the  unidirectional  flow  case.  Equation  (11-27)  with  ^rp  constant  is  the 
comparable  form  of  (11-20)  for  oscillatory  flow.  It  is  from  the  solution 
of  this  equation  that  c"(y,z,t)  is  obtained  for  determination  of  the 
longitudinal  dispersion  coefficient,  E  . 

Based  upon  the  above  discussion  it  is  concluded  that  Taylor's  analysis, 
under  the  conditions  of  uniform  flow  and  temporal  periodicity  of  u  and  c, 
is  capable  of  being  extended  for  the  determination  of  longitudinal  dispersion 
coefficients  in  oscillatory  flows  and  remains  consistent  with  the  condi- 
tions for  which  a  dispersion  coefficient  may  be  defined  as  shown  by  Aris. 
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II-D  Application  of  Dispersion  Coefficients 
Predicted  by  Taylor's  Method 

It  is  not  the  purpose  of  the  present  work  to  provide  a  detailed  anal- 
ysis of  the  conditions  under  which  values  of  E,  and  E  predicted  by 
Taylor's  method  may  be  applied  for  the  determination  of  concentration 
distributions.  However,  a  comment  regarding  this  matter  is  considered 
useful  here  in  order  to  retain  a  proper  perspective  on  the  work  presented. 

Equation  (11-12)  indicates  that  in  applying  values  of  the  dispersion 
coefficient,  steady  state  conditions  for  the  cross-sectional  mean  concen- 
tration are  not  satisfied  whereas  in  Taylor's  analysis  for  the  determina- 
tion of  the  dispersion  coefficient  it  is  assumed  that  these  steady  state 
conditions  exist.  Aris  (3)  has  shown  that  the  steady  state  conditions 
required  for  the  convective  mass  transport  to  behave  as  a  Fickian  flux  are 
approached  asymptotically  at  a  rate  determined  by  the  ratio  of  the  time 
of  convective  flow  to  the  time  of  turbulent  cross-sectional  mixing.  This 
interpretation  was  also  given  by  Taylor  (1)  on  a  qualitative  basis.  Fischer 
(21)  then  reasoned  that  a  necessary  condition  for  the  application  of  coeffi- 
cients obtained  in  this  manner  for  the  prediction  of  concentration  is  that 
c"  be  everywhere  much  less  than  c.  This  in  effect  strikes  a  compromise 
between  the  strict  steady  state  requirement  for  the  definition  of  a  dis- 
persion coefficient  and  the  conditions  for  which  it  may  be  used.  Fischer 
(21)  among  others  has  developed  detailed  criteria  for  the  application  of 
dispersion  coefficients  in  (11-12). 


CHAPTER  III 
DISPERSIVE  MASS  TRANSPORT  IN  AN  INFINITELY 
WIDE  RECTANGULAR  CHANNEL 

III-A  Steady  Unidirectional  Flow 

The  coordinate  system  has  been  selected  (Figure  2)  with  origin  at 

the  bed  of  the  channel,  z  coordinate  positive  upwards,  and  x  coordinate 

positive  in  the  direction  of  flow,  u. 


u(z) 


X z=h 


c(x, 2) 


r 


//)///// //////////^/////^/^/^^/y^y  ^ 


FIGURE  2.  DEFINITION  SKETCH  FOR  TWO 
DIMENSIONAL  SHEAR  FLOW 

The  general  procedures  used  throughout  Chapters  III  and  IV  will  be  to 

first  obtain  analytical  solutions  for  u"  and  c"  and  then  use  these  to 


calculate  the  dispersive  mass  transport  u"c"  from  which  the  dispersion 
coefficient  is  determined. 
III-A. 1   Velocity  Distribution 

It  is  assumed  that  the  turbulent  shear  stress  may  be  expressed  in 
terms  of  the  Boussinesq  approximation  and  that  the  eddy  viscosity,  e^, 
is  constant  and  equal  to  its  mean  value  over  depth.  The  equation  of 
motion  for  two  dimensional  steady  uniform  shear  flow  is  then 


p  9x    z  „,2 


(III-l) 


where  p  is  the  density  of  the  fluid  and  P  is  pressure.  It  will  be 
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assumed  in  all  analyses  that  turbulent  fluctuations  have  been  averaged 
and  incorporated  into  the  remaining  terms.  The  pressure  gradient  term 
for  this  case  may  then  be  expressed  as 

i|i=K  (III-2) 

P    oX 

where  K  is  constant.  Substituting  (III-2)  into  (III-l)  and  requiring  the 
boundary  conditions  of  no  slip  at  the  bed  and  zero  shear  stress  at  the 

surface  the  flow  boundary  value  problem  is  then 

np.  d^u  _  K  ^    ^ 

DL-  -7-—  (III-3) 

EC's:  u(0)  =  0 

dz  'z=h 
The  velocity  shear  profile  is  obtained  by  integration  of  (III-3)  and 
the  application  of  the  boundary  conditions.  The  resulting  expression  is 


K  /Z 


u(z)  =  J^(#  -  hz)  (111-4) 


■"z 
The  spatial  variation  of  the  velocity  as  defined  by  (II-2)  is  then 
obtained  by  averaging  (III-4)  over  depth  and  subtracting  out  the  depth 
mean  component  to  yield 

u"(z)  =  —  (f'  -  hz  +^')  (III-5) 

^z 

The  portion  of  the  unidirectional  flow  analysis  remaining  in  this 
section  and  portions  of  the  analyses  in  later  sections  treat  the  expres- 
sion of  the  variations  of  velocity  and  concentration  and  the  dispersion 
coefficient  as  functions  of  normalizing  parameters  other  than  the  pressure 
gradient  modulus,  K.  In  some  cases  these  alternate  formulations  are 
considered  to  be  more  useful  for  application  purposes,  whereas  in  other 
cases  they  are  presented  for  the  comparison  of  longitudinal  dispersive 
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mass  transport  in  oscillatory  and  unidirectional  flow  or  in  oscillatory 
flows  of  differing  periods  of  oscillation.  Therefore,  before  proceeding 
further  a  brief  discussion  follows  describing  these  alternate  formulations 
and  the  rationale  used  in  their  selection. 

As  noted  in  Chapter  I  the  approach  taken  by  most  investigators  in 
describing  dispersive  mass  transport  in  oscillatory  flows  has  been  to 
first  assume  a  velocity  profile  that  is  at  each  instant  in  time  similar 
in  form  to  a  corresponding  unidirectional  flow,  and  then  use  this  to 
obtain  an  expression  for  the  dispersion  coefficient,  E^,  either  directly 
through  the  method  of  moments  or  indirectly  through  Taylor's  method  of 
describing  the  convective  mass  transport  u"c"  as  a  Fickian  flux.  The 
dispersion  coefficients  thus  obtained  have  then  been  analyzed  only  as 
they  relate  to  a  "corresponding"  unidirectional  flow  coefficient.  The 
approach  taken  here  is  that  it  is  not  realistic  to  assume  a  priori  that 
the  dispersive  mechanism  in  oscillatory  flow  is  functionally  similar  to 
the  same  mechanism  in  unidirectional  floi-;  because  of  the  distinct  dif- 
ferences in  the  flow  characteristics  of  the  two  cases.  This  argument 
may  also  be  applied  to  the  more  general  case  of  comparing  dispersive 
mass  transport  characteristics  of  oscillatory  flows  of  different  periods. 
It  will  later  be  shown  that  such  comparisons  are  easily  misleading  and  if 
done  at  all  they  must  be  done  with  a  full  understanding  of  the  mechanism 
Involved.  To  illustrate  this  point,  there  are  many  ways  in  which 
unidirectional  and  oscillatory,  or  two  oscillatory  mass  transport  cases 
of  different  periods,  can  be  considered  to  be  "corresponding."  Normal- 
izing parameters  considered  in  the  present  analysis  for  comparison  of 
"corresponding"  dispersion  coefficients  are: 
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(1)  pressure  gradient  modulus,  K, 

(2)  maximum  cross-sectional  velocity,  u„,^,  which  for  the  most 
'  '  max 

general  case  of  three  dimensional  oscillatory  shear  flow  would 
represent  the  amplitude  of  the  periodic  velocity  function  located 
on  the  centerline  of  the  channel  at  the  surface, 

(3)  surface  excursion  length,  L,  experienced  by  a  particle  on  the 
free  surface  over  a  period  of  time  equal  to  one-half  the  period  of 
oscillation. 

For  the  case  of  a  unidirectional  flow,  such  as  that  being  considered 
in  this  section,  the  concept  of  an  excursion  length  as  described  above 
has  little  physical  meaning.  Thus,  excursion  length  normalization  of 
u",  c",  and  E  solutions  will  be  made  for  oscillatory  flows  only.  The 
expression  of  (III-5)  in  terms  of  u^^^  is  easily  accomplished  by 

solving  (111-4)  for  K  at  z  =  h,  or 

2e,u 
K  =  _  -^^inaii  (III-6) 

h 

which  upon  substitution  into  (III-5)  yields 

2U         2  U2 

..ii^,\  -    max  iz"-     .  ,  ,  h^x 

"  ^^^ ;T"(2-^^''3)  (III-7) 


III-A.2   Concentration  Distribution 

A  solution  for  c"(z)  as  a  function  of  K  is  obtained  by  using  the 

solution  for  the  spatial  variation  of  velocity,  u"(z),  given  by  (III-5), 

to  force  the  one  dimensional  form  of  (11-20).  Thus,  the  concentration 

variation  boundary  value  problem  may  be  stated  as 

d^c"  _  . .„  9c 


dc^ 
dz 


BC's:  z  =  0,h  ;  ^=  0 
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where, 

^  =  constant 

OS. 

u"(z)  is  given  by  (III-5). 
Substituting  for  u",  integrating  once  with  respect  to  z,  and  applying  the 
boundary  condition  at  z  =  0  produces 

dc:  ,  _l^  .  ( ll)  ^z^  _  hzi  ^  h!z^       (Ill.g) 


It  should  be  noted  that  (III-9)  implicitly  satisfies  the  boundary  condi- 

^h 
tion  at  z  =  h  since  by  definition  )      u"(z)  dz  =  0.  Equation  (III-9)  is 

then  integrated  once  more  over  z  and  the  constant  of  integration  is 

evaluated  using  the  requirement  that  c"(z)  must  have  a  zero  mean  over 

depth.  The  resulting  solution  for  the  concentration  variation  is  then 

c"(z)=^(|f)C|;-^^^-5^:    (in-10) 

^z^z 
Using  (III-6)  to  substitute  for  K  in  (III-IO)  the  expression  for 
c"(z)  as  a  function  of  u^^^^^  is 

ni\         ?1W  /9cx  pz-      hz%hfzi_h^-n      nil-ll) 

K  h^ 

z 

III-A.3   Dispersive  Mass  Transport 

The  dispersive  mass  transport  over  the  depth  of  flow  is  given  by 

m  =   \   u"c"dz  (III-12) 

0 
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Substituting  for  u"  and  c"  using  (III-5)  and  (III-IO)  respectively, 

(III-12)  becomes 

h 

.    K^  ,3c.   \  (Z'       hz  +  ^'u^  hlL^hllL.hl)  dz 
fn  =  r^%)   \  (2""''^   3^^24-  6  "^  6    45^  ''^ 

^  ^  0 

Carrying  out  the  integration  the  mass  transport  is  then 

^  =  (.  §)  .  ^]<!hl_  (III-13) 

3^    945i^R^ 

It  shall  be  assumed  for  this  analysis  and  for  the  analyses  that 
follow  that  Reynolds  analogy  relating  the  transfer  of  mass  and  momentum 
by  turbulent  processes  is  applicable  so  that 

I    =  i<  (111-14) 

z   z 

This  assumption  is  not  necessary  to  obtain  solutions  for  E|_  and  E^  using 
the  method  presented  here;  however  as  will  be  seen  later  its  use  simpli- 
fies the  forms  of  solutions  considerably  and  facilitates  interpretation 
of  the  physical  processes  involved.  Thus,  with  the  aid  of  (III-14)  and 
the  introduction  of  the  vertical  mixing  time  defined  by 

T  =hL  (III-15) 

cz 

K 

z 

Equation  (1 11-13)  may  be  written  as 

^-     2Kn3  h 
1^  =  (_  3£)       ^  cz  (III-16) 

m       V     3^;       945 

Assuming  next  that  the  required  conditions  for  defining  a  longitudinal 
dispersion  coefficient  as  discussed  in  II-B  are  satisfied,  then 

E  =  0] 

^   (-  ^)h 

or. 
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E  J-^  (III-17) 

^L    945 

Using  (III-6)  and  (III-15)  in  (III-17)  the  longitudinal  dispersion 
coefficient  as  a  function  of  u^^^^  is 

8  u^  T 

The  functional  form  of  E^  given  by  (1 11-18)  is  the  same  as  that 
obtained  by  Taylor  for  laminar  unidirectional  flow  in  a  tube, (1-8).  In 
both  cases  the  longitudinal  dispersion  varies  directly  as  the  product 
of  the  square  of  the  velocity  and  the  cross-sectional  mixing  time. 
Thus,  for  a  given  unidirectional  velocity  profile  the  longitudinal  mass 
transport  of  a  substance  increases  linearly  with  the  time  required  to 
mix  that  substance  over  the  flow  cross  section  by  diffusive  processes. 
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III-B    Oscillatory  Flow 
III-B.l      Velocity  Distribution 

The  same  assumptions  used  in  III-A  regarding  selection  of  coordinate 
system,  and  the  use  of  the  Boussinesq  approximation  for  expressing  the 
viscous  term  in  the  equation  of  motion  will  be  applied  here.     The  govern- 
ing equation  of  fluid  motion  for  uniform  unsteady  two  dimensional   shear 
flow  is  then 

3u  _      1  9P   .  p    _3!u  (III-19) 

3t  "  ■  P  "^^T        z     8z2 

For  the  oscillatory  flow  case  it  will  be  assumed  that  long  wave  phenomena 

are  of  primary  interest  and  therefore  the  pressure  distribution  over  depth 

is  hydrostatic.  It  is  then  reasonable  to  assume  that  the  pressure  gradient 

in  (II 1-19)  is  periodic  and  of  the  form 

i|P  =  Ke^'^^  (III-20) 

p  8x 

where  K  is  constant,  as  before,  and  a  is  the  angular  frequency  of  oscilla- 
tion. Note  that  only  the  real  part  of  (III-20)  has  physical  meaning. 
To  obtain  a  solution  assume  that  u(z,t)  is  periodic  in  time  and  of 

the  form 

u(z.t)  =  U3(z)e''"'^*  (111-21) 

where  u  (z)  is  a  complex  function.  The  boundary  value  problem  is  then 

formulated  by  incorporating  (III-20)  and  (III-21)  into  (III-19)  which 

along  with  the  boundary  conditions  of  no  slip  at  the  channel  bed  and 

zero  shear  stress  at  the  surface  becomes 

d^u 
DE:  E  -^-  iau,  =  K  (III-22) 

—   z  dz^      s 

EC's:  u^(0)  =  0 


du 


dz  iz=h 


=  0 
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This  problem  has  been  previously  solved  by  Lamb  (18)  and  Segall  (19) 

who  obtained 

t     *^   iKn   cosh  g(l+i)(h-z)nj'at       (III-23) 
u(z,t)  =  —  LI  -  cosh  3(l+i)h  ^^ 

where ,  

with  only  the  real  part  of  (III-23)  having  any  physical  meaning.  Its 
form  is  clearly  that  of  a  damped  progressive  shear  wave  propagating 
upward  through  the  water  column  and  causing  a  temporal  phase  shift  in 
the  velocity  as  a  function  of  z.  It  was  also  pointed  out  by  Segall  and 
Gidlund  (20)  that  the  solution  as  expressed  by  (III-23)  correctly  predicts 
flow  reversals  in  the  lower  momentum  layers  of  the  fluid  near  the  bed 
prior  to  a  shift  in  flow  direction  of  the  higher  momentum  layers  near 
the  free  surface.  These  effects  are  extremely  important  to  the  disper- 
sive transport  mechanism. 

The  spatial  variation  of  the  velocity,  u"(z,t),  is  obtained  by 
averaging  (III-23)  over  depth  and  subtracting  out  the  mean  from  the 
total  velocity  according  to 

u"(z,t)  =  u(z,t)  -  u(t) 

The  result  is  then 

u"(z.t)  =  -^^ {  ^^'"^^^^"^^^  -  cosh  3(l+i)(h-z)}   (in-25) 

a   cosh  6(l+i)h    B(l+i)h 

For  reasons  that  will  become  obvious  in  the  next  section  the  velocity 
variation  solution  as  given  by  (III-25)  will  now  be  expanded  in  a  Fourier 
cosine  series  of  the  form 

oo 

u"(z,t)  =  y  a„  cosl^e^'^*  (III-26) 


n=l 
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where  coefficients  a^^  are  complex.  The  a^  term  in  the  series  has  been 
omitted  to  satisfy  the  zero  mean  requirement  for  u".  It  is  noted  by 
Hildebrand  (22)  that  an^  piecewise  differentiable  function  may  be  com- 
pletely represented  by  a  cosine  series  of  the  form  a^  cos  -^    over  the 
closed  interval  0  <  z  <  h.  Proceeding  then, 

2   r      iK  sinh  B(l+i)h      mrz  , 
\  =  h       \        a3(l+i)h  cosh  3(1+1  )h  "'  h 

^0 

The  first  term  on  the  right  hand  side  integrates  to  zero  while  the  second 
term  must  be  integrated  by  parts  twice  which  yields  for  a^, 

.     -      ^K  3h(l-i)  sinh  g(l+i)h fTTT  ?7^ 

\  -   a  cosh  3(l+i)h  L(n7T)2  +  iZ(3h)-J         (III-27) 

Since  only  the  real  part  of  u"(z,t)  is  of  any  interest,  (III-27)  will  now 
be  put  in  polar  form  for  incorporation  into  (II 1-26).  Considering  each 
complex  term  separately, 

(1-i)  =  /Z    e'^l  (III-28) 

sinh  3(l+i)h  =  ^  (cosh  23h  -  cos  23h)''^  e^"^'^        (III-29) 

cosh  3(l+i)h  =  "^  (cosh  23h  +  cos  23h)''^  e"*^^        (III-30) 

C(nTr)2  +  i2(3h)^]  =  [(nTr)"*  +  4(3h)'*f^  e'"'"^  (III-31) 

where, 

-  +=.n-^  ,cosh3hsin_3h\  fTII-32^ 

^'   =  ^^"   ^sinh  3h  cos  3h^  ^^"  ^^^ 

4.,«-^  /Sinh  Bh  sin  ghx  fIII-33) 

«2  =  t^"   (cosh  3h  cos  3h^  ^     ' 

0,3  =  tan-^  L^^^^-}  (ni-34) 

(nvr)^ 
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To  express  u"(z,t)  in  its  final  form  it  is  again  assumed  that 
Reynolds  analogy  is  applicable  allowing  the  use  of  (III-15).  A  non- 
dimensional  time  relating  the  period  of  oscillation  to  vertical  mixing 
time  is  then  introduced  by 

T«  =  _£Z  (III-35) 

'z   T 

With  the  aid  of  Reynolds  analogy  and  (III-15),  (III-24),  and  (III-35) 

the  3h  arguments  appearing  throughout  a^  become 


3h  =Vti1' 


(III-36) 


The  final  form  for  u"(z,t)  as  a  function  of  the  pressure  gradient  is 
arrived  at  by  combining  (III-26)  through  (III-36)  to  yield. 


2T'    (cosh  2VVr   -  cos  2VuT')     j, 
u"(z,t)  =  KT  {-^ z=^ 7=^  } ' 


00 

I 

n=l 


IT  (cosh  2v^TrT^  +  cos  2  v  ttT^ 

r-r..  fi^z  ^i(at  -  V  +  "1  -  o'z     -  ct3(n)) 
cos  -fp  e 


[(nTT)-*  +  (2ttT;)2J 


(III-37) 


where , 


ai  =  tan        { 


cosh 


ttT'   sin  V  ttT' 


sinh    a/ttTI  cos  v/ttT^ 


a2  =  tan         {- 


sinh 


''ttT^  sin  a/tttZ 


cosh 


-1    2TrT' 
a3(n)  =  tan     (r^z) 


ttT'   cos  ^/uP 


(III-38) 
(III-39) 
(III-40) 


Expressions  for  u"(z,t)  as  a  function  of  the  surface  amplitude, 

u  of  the  oeriodic  velocity  function  and  the  surface  excursion  length, 

max'  ^ 

L,  experienced  by  a  particle  during  one-half  of  the  period  of  oscilla- 
tion are  obtained  by  solving  for  K  as  a  function  of  u^^^  and  L  using 
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(III-23).  Details  of  this  development  are  presented  in  Appendix  A,  the 
results  of  which  are 


\P[    TT{cosh  2  v^ttT;  +  cos  2  ^A(V) ' 


—  u 


T  (cosh  V  ttT^  -  cos 


Vr^) 


max 


(111-41) 


and 


K  = 


v^  TT^   L  (cosh  2  v^  +  cos  2  aAtTP' 


(III-42) 


T^  (cosh  -/ttT;  -  cos 


.t;; 


Substitution  for  K  in  (III-37)  using  these  two  expressions  then  yields 

2u   {ttT'  (cosh  2  y^  -  cos  2  v^)}^^ 
u"(z,t)  =  - 


max 


00 

I 

n=l 


(cosh   a/Tt^  -  cos   VttT^) 
n-TTZ  „i(at  -  T-  +  ai  -  ct2  -  013) 


cos-pe 


;iII-43) 


C(nTr)^  +  (2ttT;)^]^ 


and. 


u"(z,t)  =  - 


2TrL 


{  uT'    (cosh  2v/7r!  -  cos  2  ^At!)}^" 


T  (cosh  vmtTI  -  cos  VttT' ) 


nTTZ  „i(at  -  4-  +  ai  -a2  -  aa, 
cos  -r—  e   '  ^  z  d, 


n=l 


\_{mY  +  (2ttT;)^T 


(III-44) 


III-B.2      Concentration  Distribution 

For  temporally  periodic  u  and  c,  uniform  flow,  and  constant  eddy 
diffusivity  it  was  shown  in  II-C  that  the  transport  diffusion  equation 
as  described  from  a  coordinate  system  moving  with  the  cross-sectional 
mean  velocity  could  be  written  as 


9C         |,     9  c       _       II   8c 
3t    "    ^z  9Z2       '  ■"     3^ 


(III-45) 


46 


where  the  rate  of  change  with  respect  to  5  will  hereafter  be  implied. 
The  concentration  variation,  c"(z,t),  shall  be  assumed  to  be  of  the  form 

00 

c"{z.t)  =  c^(z)e^"^*  =  ^a;cos2-e^--t  (iii-46) 

m=l 

where  a'  is  complex.  This  expression  for  c"  is  seen  to  satisfy  the  require- 
m 

ments  of  periodicity  and  zero  diffusive  flux  of  substance  across  the  flow 

boundaries.  Moreover,  if  (111-46)  is  substituted  into  (III-45)  and 

operated  on,  the  left  hand  side  is  composed  of  cos  ^  terms  only.  Thus, 

to  satisfy  the  conditions  of  equality,  the  expansion  of  (II 1-25)  for  u" 

in  a  cos  ^series  is  justified.  The  formal  statement  of  the  boundary 
h 

value  problem  for  c"(z)  is  then 

where  |t  is  constant  and  u"(2)  is  defined  by  (III-26).  Substitution  for 

c"  and  u"  transforms  (II 1-47)  to 
s     s 


oo 

I 


iK^¥  ^  '^K  "^  T  =  <-  If)  •  Y.\  "^  T 


m=l  n=l 

which  due  to  the  orthogonal  properties  of  cos  ^  and  cos  -p-  requires 
that  m=n  for  a  non-trivial  solution.  The  Fourier  coefficient  a^  can 
then  be  written  as 

a-  =  (_  9£)  ^ , (III-48) 

"^  ^  '^^  R  r(rM)^  .  i(g^^)] 

^z 

The  required  solution  for  c"(z,t)  is  arrived  at  by:  (1)  using 
(III-26)  and  (III-37)  to  determine  the  polar  form  of  a^^;  (2)  transforming 
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the  denominator  of  (II 1-48)  to  polar  form;  (3)  applying  the  definitions 

for  T      and  T'  where  necessary;  and  (4)  substituting  the  resulting  expression 
cz  z 

into  (III-46).     The  variation  component  of  concentration  as  a  function  of 

the  pressure  gradient  is  then, 

3     2(cosh  2  ^/Tr   -  cos  2  v^^)     ,. 


3c. 


c"(z.t)  =  (-  f|)  KTn;5  { 


} 


m=l 


TT(cosh  2'v/^  +  cos  2  v^ttT^) 


muz  „i(at  -  ^  +  ai  -  a2  -  2a3{m)) 
cos  -T—  e  ■* 

C(mir)'*  +  (2uT;)2] 


(III-49) 


where  ai,  a^,  and  aaCm)  are  given  by  (III-38)  through  (III-40)  with  m 
replacing  n  in  (III-40). 

Solutions  for  c"(z,t)  as  a  function  of  u^^^,  and  L  are  obtained  by 
applying  (III-41)  and  (III-42)  to  (III-49),  to  yield 

(cosh  a/tttJ  -  cos    vVrJ) 

mrrz  „i(at  -  ^  +  ax   -  ctz   -  2a3(m)) 
cos  — r—  e  ■♦ 

h 


[(rnir)'^  +  (20^] 


(III-50) 


m=l 


and 


c"(z,t)   =  2L(7rTp2    (||) 


(cosh  2  v^  -  cos  2  y^)^' 
(cosh  v/uT^  -  cos  v^^) 


^^3  mjz  gi(at--+a, 
n  


aa 


2a3(m) 


m=l 


C(miT)'*  +  (2uT;)^] 


(III-51) 


III-B.3   Dispersive  Mass  Transport 

For  the  oscillatory  flow  case  the  mass  transport  considered  will  be 
the  time  mean  transport  averaged  over  one  period  of  oscillation.  Thus, 


48 


:ni>^=A  R,{u")   •R3(c")dzV 


(III-52) 


the  dispersive  mass  transport  is  given  by 
<m>. 

;o 

Using  (III-37)  and  (III-49)  for  u"(z,t)  and  c"(z,t),  the  integra- 
tion over  depth  of  the  product  Rg(u")  •  Rg(c")  reduces  the  implied 
double  summation  over  m  and  n  to  a  single  sum  through  the  orthogonal 
properties  of  the  functions  cos  ^    and  cos  ^    since, 


^«„    HTTZ    ^„^    mTTZ 

COS  -fp  cos  -^ 


dz  = 


^  ,   n  =  m 


'0  to  >  n  ?*  f" 

Performing  this  integration  and  averaging  the  results  over  T  yields. 


<m>. 


^-  h  Knn'2       (cosh  2  yTr;  -  cos  2  a/^) 

dC\ Z       ,    Z i. 

■b^'       27r 

cos  a3(n) 


(cosh  2  vttTT  +  cos  2  vttT^) 


^z^'CCnTT)'*  +  (27rr)2]2 


(III-53) 


The  longitudinal  dispersion  coefficient  is  next  introduced  using 
Taylor's  analogy  defining  the  convective  mass  transport  <m>j  as  a 
Fickian  flux  so  that 


<m>,     =  <-  E^  h(||)>, 


(III-54) 


Noting  from  (III-40)  that 

cos  aa   (n)  = 


(nu)' 


[(n^)"*  +  [i^vvy 
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and  combining  (III-53)  and  (III-54)  the  time  mean  longitudinal  dispersion 
coefficient  as  a  function  of  the  pressure  gradient  is  expressed  as 


<c  >  =  K^m^^   (cosh  2  V^  -  cos  2  a/Tt^) 


^x^T 


TfT 


n=l 


(cosh  2  vVi^  +  cos  2  v^rT^) 


(htt)' 


.  C(n7r)'*+(2TTT')2] 


(III-55) 


Solutions  for  <E  >-r  as  a  function  of  u_,^  and  L  may  be  obtained  by 

X  I  illdA 

substituting  directly  for  K^  in  (III-55)  using  (III-41)  and  (III-42). 
The  resulting  expressions  are  as  follows: 

(cosh  2  VTr  -  cos  2  V^) 


'h>l 


■^uL  T  T- 
max   z 


(cosh 


(nTT)  = 


ttT'  -  cos 

z 


vr;)^ 


^^1^  L{mrH2i^V^)'T 


(III-56) 


and. 


TT^  L^T'"  (cosh  2y7r  -  cos  2  ^Aiv) 


X  T 


(cosh 


ttT'  -  cos  v^^)^ 


iiml: 


[(mr)^  +  (27Tr  )^T 


(III-57) 


Equation  (III-56)  is  considered  to  be  the  most  useful  and  descrip- 
tive of  the  solutions  presented  for  <E^>j  in  two  dimensional  oscillatory 
shear  flow.  Therefore,  to  facilitate  a  discussion  of  these  results  in 
the  next  section,  (1 11-56)  is  non-dimensional ized  by  defining 

<E  > 


^;  = 


-X  T 


max 


(III-58) 
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and  restating  it  as 


E".  = 


TT  T^Mcosh  2    V  ttT'     -  cos  2    'v/  ttT'] 


(cosh    v^ttTT     -  cos    \/  ttT  ' )  2 


00 


(n^)' 


^j^Y"  [(n^)'  +  (2^t;)2]^ 


(111-59) 
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III-C    Discussion  of  Two  Dimensional  Shear  Flow  Results 

III-C.l      Dispersion  in  Unidirectional   Flow  as  a  Limit  of  the  Oscillatory 
Flow  Case  ~~~ 

This  discussion  shall  begin  by  relating  the  oscillatory  flow  disper- 
sion coefficient,  <E  >-.,  to  its  "corresponding"  unidirectional   flow 
coefficient,  E,  ,  as  has  been  previously  done  by  Holley  et  al . ,  Awaya, 
and  Fukuoka;  and  shall   then  proceed  to  show  how  the  unidirectional   process 
is  in  fact  a  limiting  case  of  an  oscillatory  process  which  has  its  own 
distinct  characteristics.     Thus,  solutions  obtained  in  III-B.3  for 

<E  >-r  as  a  function  of  the  pressure  gradient  and  the  surface  amplitude 
X  T 

of  the  velocity  are  normalized  by  the  "corresponding"  solutions  for  E^^ 
obtained  in  III-A.3.     Dividing  (III-55)  by  (III-17)  and  (III-56)  by 
(III-18)  yields 


<E^>j      g45       (cosh  2    V^ttT^  -  cos  2 


^L        ^^"""z     (cosh  2    vXTr;  +  cos  2  v^^^) 


}     Ml (III-60) 

^[(nTT)^  +  (2ttT;)^]^ 


for  the  same  K,  and 


<E  >j       9457tT^       (cosh  2    ^/\l^  -  cos  2  VttT  ) 
^L  ^  (cosh    v/^  -  cos    v^^)2 


00 

z 


(nu)' 


(III-61) 


^    [{mr  +  (2ttT;)^]^ 

for  the  same  u   .  All  calculations  involving  infinite  summations  for 
max 

two  dimensional  shear  flow  dispersion  use  200  terms  in  each  sum,  well 
within  fourth  place  accuracy. 
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The  expressions  given  by  (III-60)  and  (III-61)  are  plotted  in  Figure 
3  vs.  the  non-dimensional  time  1/T^  corresponding  to  T'  defined  by 
Hoi  ley  et  al.  Also  included  in  Figure  3  is  the  solution  obtained  by 
Hoi ley  et  al.  for  <E  >^/E,  for  the  oscillatory  linear  velocity  profile 
discussed  in  I-B.2.  It  is  noted  that  all  three  solutions  approach  the 
limiting  value  of  one-half  for  situations  when  the  period  of  oscillation 
is  greater  than  the  vertical  mixing  time  in  accordance  with  the  findings 
of  Okubo  and  Bowden.  However,  the  difference  in  behavior  of  the  two 
solutions  given  by  (III-60)  and  (III-61)  is  dramatic  and  immediately 
raises  the  question  of  whether  or  not  a  plot  such  as  Figure  3  is  the 
most  meaningful  method  of  illustrating  the  characteristics  of  dispersive 
mass  transport  in  oscillatory  flow.  For  the  case  of  the  pressure 
gradient  normalized  solutions  the  ratio  ^E/j/^l  decreases  rapidly 
with  increasing  T'  because: 

(1)  For  unidirectional  flow  the  pressure  gradient  and  viscous 
forces  are  in  equilibrium  so  as  to  produce  a  constant  unidirec- 
tional shear  flow  for  the  dispersive  transport  of  substance. 

(2)  In  oscillatory  flow  the  pressure  gradient  is  assumed  to  be 
simple  harmonic  in  time  and  therefore  is  in  constant  balance 
with  the  varying  inertia  and  friction  forces.  As  the  period 

of  oscillation  decreases  the  inertial  effects  become  very  large 
such  that  in  the  limiting  case  little  or  no  flow  would  be 
Induced.  This  results  in  little  or  no  dispersive  transport.  This 
effect  of  the  period  of  oscillation  on  velocity  can  be  seen 
from  (II 1-23). 
The  solutions  normalized  by  u„,,  exhibit  a  much  slower  decrease  of 

maX 

<E  >-[./£.  with  increasing  T^  simply  because  in  this  case  the  surface 
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FIGURE  3  NON-DIMENSIONAL  DISPERSION  COEFFICIENT 
<E  >j/E|_,  AS  A  FUNCTION  OF  T' 
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velocities  of  the  unidirectional  and  oscillatory  flow  are  required  to  be 
the  same.  In  terms  of  pressure  gradients  this  would  correspond  to  the 
situation  where 

osc    umdir 
Based  upon  these  results  it  becomes  apparent  that  due  to  the  very 

different  hydrodynamic  characteristics  of  unidirectional  and  oscillatory 
flow  systems,  the  interpretation  of  the  nature  of  dispersive  mass  trans- 
port in  an  oscillatory  system  as  compared  to  the  same  process  in  a  uni- 
directional system  is  misleading  and  can  vary  widely  depending  upon  how 
one  chooses  to  relate  the  two  systems.  This  is  demonstrated  by  Figure  4. 
In  this  figure  the  longitudinal  dispersion  coefficient,  <Ex^j.  as  given 
by  (III-56)  is  plotted  against  the  vertical  mixing  time,  T^^,  for  four 
periods  of  oscillation  ranging  from  22,536  to  223,560  seconds  and 

a  u   of  1  ft/sec.  Also  plotted  is  the  solution  for  the  unidirectional 
max 

flow  coefficient,  E,_,  as  given  by  (III-18).  As  shown  by  Figure  4  the 
behavior  of  E,  and  <E  >j   are  clearly  different.  The  unidirectional  flow 
coefficient  varies  directly  with  the  vertical  mixing  time  which  for  little 
or  no  turbulent  mixing  over  the  water  column  allows  the  shear  flow  to 
transport  higher  concentrations  of  substance  far  downstream.  The  behavior 
of  the  oscillatory  flow  coefficient,  however,  is  governed  by  a  type  of 
resonant  interaction  between  the  period  of  oscillation  and  the  vertical 
mixing  time.  As  the  period  of  oscillation  is  increased  the  resonant 
peak  shifts  to  the  right  and  likewise  increases  until  the  limiting  case 
is  reached  where  the  peak  is  infinitely  large  and  the  values  predicted 
for  <E  >j   are  exactly  one-half  those  predicted  for  the  unidirectional 
flow  coefficient  E|^. 
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The  resonant  characteristic  of  <£^>j  may  also  be  used  to  explain 

the  abrupt  decrease  in  the  ratio  <£y>j/\   as  shown  in  Figure  3.  The 

value  of  T'  at  which  the  knee  of  the  curve  occurs  in  Figure  3  corresponds 

to  the  value  of  T   in  Figure  4  at  which  <E  >t  for  a  given  T  begins  to 
cz  «  ' 

diverge  from  the  unidirectional  solution.  At  this  point  E|_  continues 

to  increase  whereas  <E  >^  peaks  and  begins  to  decrease. 

III-C.2   Characteristics  of  Dispersive  Mass  Transport  in  Oscillatory  Flows 

The  physical  reasoning  behind  the  resonant  behavior  of  <E^>j  is 
straightforward  and  is  begun  by  considering  the  extreme  cases  of  T  «  T^^ 

and  T  »  T^^: 

(1)  T  «T      -   In  this  case  the  degree  of  vertical  mixing  is  so 
small  during  one  period  of  oscillation  that  an  elemental   volume  initially 
residing  in  the  water  column  at  elevation  z,  and  containing  an  initial 
concentration  c,  would  remain  at  this  elevation  thus  being  transported 
over  the  closed  pathline  of  flow,  and  returned  to  its  initial   position 
with  no  net  longitudinal   dispersion  having  occurred.     Conversely.in  a 
unidirectional   flow  with  little  or  no  vertical  mixing  the  longitudinal 
dispersion  would  be  yery  large. 

(2)  T  »  T      -  Here  the  rate  of  vertical  mixing  is  so  rapid  that 


cz 


there  is  no  time  for  the  velocity  shear  profile  to  transport  the  sub- 
stance longitudinally  before  it  loses  its  identity  through  vertical 
mixing.  In  this  case  the  oscillatory  and  unidirectional  flow  dispersive 
processes  behave  in  a  similar  manner  and  are  both  small. 

Thus,  it  is  seen  that  for  an  oscillatory  flow  the  longitudinal 
dispersive  mass  transport  becomes  small  for  both  T  «  T^^  and  T  »  T^^, 
whereas  for  a  unidirectional  flow  the  longitudinal  dispersive  mass  trans- 
port varies  directly  with  T^^.  This  resonant  nature  of  <E^>-p  can  be 


57 


further  explained  on  physical  grounds  by  beginning  with  Case  1  above  and 
gradually  allowing  the  vertical  mixing  time  to  decrease.  In  this  example 
the  velocity  shear  profile  in  conjunction  with  the  increased  vertical 
mixing  causes  a  net  longitudinal  mass  transport  to  occur  over  a  period 
of  oscillation.  This  effect  continues  to  increase  as  T^^  is  decreased 
until  the  optimum  ratio  between  the  vertical  mixing  time  and  the  period  of 
oscillation  is  reached.  At  this  point  the  longitudinal  dispersive  mass 
transport  has  reached  its  maximum  value.  A  further  decrease  in  T^^ 
begins  to  introduce  the  effect  noted  in  Case  2  where  the  vertical 
mixing  is  now  too  rapid  thus  causing  <^^>j   to  decrease. 

Additional  physical  insight  into  the  behavior  of  <E^>j  can  be  ob- 
tained by  considering  the  relative  longitudinal  dispersion  associated 
with  two  oscillatory  flows,  one  having  twice  the  period  of  oscillation 
of  the  other.  The  two  periods  considered  are  44,712  and  89,424  seconds 
corresponding  to  the  periods  of  a  semi-diurnal  and  diurnal  tide 
respectively. 

The  first  of  the  two  comparisons  is  illustrated  by  Figures  5a  and 

5b.  Once  again  the  expression  used  for  <^^>j   is  given  by  (III-56)  with 

u   =1  ft/sec  for  both  the  semi-diurnal  and  diurnal  cases, 
max 

As  shown  in  Figure  5a  after  an  elapsed  time  of  t  =  89,424  sec  the 
semi-diurnal  tide  will  have  completed  two  full  cycles  of  excursion 
length  L  ,  while  the  diurnal  tide  will  have  completed  one  full  cycle  of 
excursion  length  L^  =  2L^.  In  Figure  5b  the  favorable  effect  of  the 
longer  excursion  length  on  longitudinal  dispersion  is  seen  for  flows  of 
different  periods  of  oscillation  with  the  same  u^^^.  The  second  compari- 
son betv/een  the  semi-diurnal  and  diurnal  period  flows  uses  excursion 
length  as  a  normalizing  parameter  and  is  illustrated  by  Figures  6a  and 
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6b.     In  this  comparison  the  velocity  for  the  diurnal  tide  is  scaled 
downward  so  as  to  produce  the  same  excursion  length  as  the  semi-diurnal 
tide.     Figure  6a  shows  that  after  an  elapsed  time  of  t  =  89,424  sec  the 
semi-diurnal  tide  will  have  completed  two  full  cycles  of  excursion 
length  L  as  before.     However,  the  diurnal  tide,  because  of  the  reduced 
velocity,  has  completed  one  full  cycle  of  the  same  excursion  length. 
Values  of  <E  >-.  are  plotted  vs.  T      using  (III-57)  in  Figure  6b.     As 
seen  from  this  plot  the  semi-diurnal  dispersion  has  remained  unchanged 
from  Figure  5b  whereas  the  diurnal  curve  is  now  everywhere  less  than 
the  semi-diurnal  curve.     This  is  not  surprising  since  from  (III-57)  the 
fixing  of  excursion  length  and  vertical  mixing  time  leaves  only  the 
decrease  in  the  period  of  oscillation  as  a  means  of  increasing  longitudinal 
dispersion. 

The  results  of  this  analysis  on  the  predicted  behavior  of  <E^>-|. 
in  two  dimensional  oscillatory  shear  flow  are  summarized  in  non- 
dimensional  form  by  Figure  7.     In  this  figure  the  non-dimensional 
longitudinal  dispersion  coefficient  as  defined  by  (III-59)   is  plotted 
vs.  T'.     Thus,  by  knowing  the  period  of  oscillation,  the  surface 
amplitude  of  velocity,  the  water  depth,  and  the  vertical  eddy  diffusi- 
vity,  one  could  use  Figure  7  to  obtain  a  predicted  value  for  the 
longitudinal  dispersion  coefficient.     The  resonant  value  of  T^  for 
which  <E  >j  is  a  maximum  is  seen  from  Figure  7  to  be 

T-   =  1.58  (III-62) 

This  relationship  is  applicable  to  all  oscillatory  flows  in  in- 
finitely wide  rectangular  channels  provided  the  assumptions  made  in  this 
analysis  are  reasonably  well   satisfied. 

A  look  at  the  kinematic  structure  of  the  flow  field  and  the 
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FIGURE  6  COMPARISON  OF  DISPERSION  FOR  DIURNAL  AND  SEMI-DIURNAL 
PERIODS  OF  OSCILLATION  WITH  SAME  SURFACE  EXCURSION 
LENGTH 
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FIGURE  7  NON-DIMENSIONAL  LONGITUDINAL  DISPERSION  COEFFICIENT 
E',  AS  A  FUNCTION  OF  T^  FOR  TWO  DIMENSIONAL 

OSCILLATORY  SHEAR  FLOW 
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concentration  distribution  over  depth  is  presented  for  three  values  of 
T'  in  Figures  8,  9,  and  10.  In  these  figures,  non-dimensional  forms  of 
the  solutions  for  u"(z,t)  and  c"(z.t)  given  by  (III-43)  and  (III-50) 
respectively  are  plotted  over  the  non-dimensional  depth,  z/h.  The  non- 
dimensional  forms  used  are 


ND   u_„ 
max 


c 


(III-63) 


(III-64) 


8^^  max 

Each  figure  shows  the  u"  and  c"  profiles  for  maximum  flow  and  slack 
water.  Several  interesting  points  are  illustrated  here.  First,  the 
superior  correlation  of  u"  and  c"  over  depth  in  Figure  9  when  compared 
to  the  same  correlation  in  Figures  8  and  10  demonstrates  qualitatively 
why  <E  >j   is  greater  for  T^  =  1.58  than  for  the  other  values  of  T^ 
considered.  Second,  the  phase  dependency  of  the  velocity  with  depth  is 
clearly  shown  in  all  figures  near  slack  water  with  flow  reversals 
occurring  in  the  lower  momentum  layers  near  the  channel  bed  prior  to 
a  shift  in  the  upper  higher  momentum  layers.  Finally,  it  is  interesting 
to  note  that  the  shear  profile  is  confined  to  a  decreasingly  thinner 
layer  near  the  bed  with  the  upper  portion  of  the  profile  becoming 
flatter  as  T'  is  increased  (increasing  eddy  viscosity  by  Reynolds 
analogy). 
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CHAPTER  IV 
DISPERSIVE  MASS  TRANSPORT  IN  RECTANGULAR 
CHANNELS  OF  FINITE  WIDTH 

IV-A  Oscillatory  Flow 

The  physical  problem  and  the  general  analytical  approach  used  to 

solve  it  are  the  same  in  this  chapter  as  those  discussed  in  Chapter  III 

with  the  added  feature  of  a  finite  channel  width.  The  effect  of  this 

feature  is  to  introduce  a  three  dimensional  shear  flow  with  vertical  and 

transverse  velocity  variations  which  are  used  to  force  the  transport 

diffusion  equation  to  obtain  a  solution  for  c"(y,z,t).   Solutions  for 

u"(y,z,t)  and  C'ly.z.t)  are  then  correlated  over  the  flow  cross  section 

and,  in  the  case  of  oscillatory  flow,  averaged  over  the  period  of 

oscillation  to  obtain  the  longitudinal  dispersive  mass  transport. 

The  coordinate  system  selected  has  its  origin  at  the  center  of  a 

rectangular  section  whose  sides  correspond  to  z  =  ±h,  and  y  =  ±w/2.  As 

shown  in  Figure  11,  the  open  rectangular  channel  is  mathematically 

represented  by  the  lower  half  of  the  section.  Selection  of  the  coordinate 

system  in  this  manner  preserves  the  syrmetry  of  the  problem  about  the 

origin  and,  as  will  be  shown  later,  correctly  predicts  the  same  dispersive 

mass  transport  for  equal  degrees  of  skewness  in  channel  geometry  in  either 

the  vertical  or  lateral  directions. 
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FIGURE  11.    DEFINITION  SKETCH  FOR  THREE 
DIMENSIONAL  SHEAR  FLOW 

IV-A.l        Velocity  Distribution 

Using  the  same  assumptions  that  were  made  in  Chapter  III  regarding 

eddy  viscosity,  time  averaging  of  turbulent  fluctuations,  and  the  form 

of  the  pressure  gradient  term,  the  equation  of  motion  for  uniform  flow 

in  the  x  direction  may  be  written  as 


|u  ._  .Keiat  ,  -     8!u  ,  -     31U 

at  z  3^2       y  9y2 


If  a  solution  of  the  form 


u(y,z,t)  =  u  (y,z)  e 


iat 


(iv-i: 


(lV-2) 


is  assumed  to  be  valid  for  (IV-1)  then  the  three  dimensional  shear  flow 
boundary  value  problem  for  oscillatory  flow  may  be  formally  stated  as: 


PDF: 


3'u. 


'^  8z^ 


3'U. 


+  e 


y  9y^ 


iau.  =  K 


{IV-3) 


BC's:  u^(y,±h)  =  0 
u^ltf,.)  -  0 

It  is  noted  that  physical  considerations  require  that  u^(y,z)  be 
an  even  function  in  y  and  z.  Thus,  a  solution  to  (IV-3)  of  the  form 
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00       CO 

m=0  n=0 
is  assumed.  An  examination  of  (IV-4)  with  the  statement  of  the  boundary 
value  problem  above  shows  that  the  assumed  solution  satisfies  the  no 
slip  conditions  on  the  perimeter  of  the  section.  Moreover,  (IV-4)  also 
satisfies  the  condition  of  zero  vertical  shear  at  the  free  surface,  and 
zero  lateral  shear  at  the  centerline  of  the  channel.  The  assumed  solution 
is  then  substituted  into  (IV-3)  and  the  Fourier  coefficients,  a^^,  are 
determined  in  the  usual  manner  through  the  use  of  the  orthogonal  proper- 
ties of  the  cosine  function.  Details  of  this  analysis  are  presented  in 
Appendix  B.  Solving  for  a   and  thus  determining  the  form  of  Ug(y,z), 
the  complete  solution  for  the  velocity  function  u(y,z,t)  using  (IV-4) 
and  (IV-2)  is  then 


00     00 


To  obtain  the  spatial  variation,  u"(y,z,t),  the  cross-sectional  mean 
velocity,  u(t),  must  be  subtracted  out  of  (IV-5).  The  mean  component  of 
u  is  defined  by  h    w/2 

J-h   ^-w/2 

Performing  the  integration  in  (IV-6)  and  subtracting  the  result  from 

(IV-5)  yields 

00       00 

u"(y.z.t)  =  y  ya„„  {COS  i^^  cosi^^  -  g(m.n))e-t    (IV-7) 
m=0     n=0 
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where , 


and. 


g(m,n)  =  ^^^ (IV-8) 

7T2(2m+l)(2n+l) 


16K  (-1)"^"        (iv-9) 


"""    uM2m.l)(2n.l)  (  I    L^^f  *  S  L^^^^f  *   i") 


Since  only  the  real  portions  of  these  expressions  have  any  physical 
meaning,  (IV-7)  is  transformed  to  polar  form  as  was  done  in  Chapter  III 
for  the  two  dimensional  case.  In  addition,  it  will  be  assumed  once  again 
for  simplification  of  results  that  Reynolds  analogy  applies  thereby 
allowing  the  introduction  of  the  vertical  mixing  time  T^^  and  the  defini- 
tion of  a  comparable  lateral  mixing  time  as 

T    =  i^:M'  (iv-io) 

The  solution  for  u"(y,z,t)  as  a  function  of  the  pressure  gradient  may 
then  be  written  as 

/■,Nm+n  gi(at-ai(m,n)) 
u"(y,z,t)   -'^    ^ 


m  =  0    n  =  0     (2  m+l){2  n+1)  Rf  (m,n) 


(2_m±l)iz      cos  ilililM  .  g(m,n)  I  (IV-U) 


cos  2h 


where , 


R, 


d<n.n)A^l^^^r*^^^^^y\       *1  ("-12' 


SIT^^       2        -■       2.T^y^       2 


cCm.n)  =  tan-'  J     — ^(2  m-Uji,-,,! T      A2  n^l)^.,     I        ""-"' 

2^^         2       ]"2.T^/       2         J 
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In  order  to  put  (IV-11)  in  a  more  usable  form,  (IV-5)  is  used  to 

obtain  a  relationship  between  the  pressure  gradient,  K,  and  u^^^,  where 

u   is  now  defined  as  the  amplitude  of  the  periodic  velocity  function  at 
max 

the  free  surface  centerline.  Thus,  from  (IV-5) 

u(o.o,t)  =^   y  . 


m  =  0    n  =  0 


■16  K(-ir^"  e"*'^^ 


TrM2  m+l)(2  n+1)   {e^[ll^]^+  ^^[11^]^+  ia} 


:iV-14) 


which  upon  maximization  with  repsect  to  time  yields  the  relationship 

.3 


K  = 


U         TT" 

max 
8T(J) 


(IV-15) 


where. 


y^    V^(-ir"(R,(m.nl-1)'-' 
:fr-o    fm  (2m+l)(2  n+l)R,(m,nJ 

00  oo 


m  =  0     n  =  0 


(2  m+l)(2  n+1)   Ri(m,n) 


(IV-16) 


The  details  of  developing  (IV-15)   from  (IV-14)  are  presented  in  Appendix 
C.     Substitution  of  (IV-15)   into  (IV-11)   for  K  yields  the  following 
expression  for   u"  as  a  function  of  u^^^ 


u"(y,z,t) 


max 


/   -.xm+n  i(at-ai(m,n)) 
*        fr^    ^T-Q     (2m+l)(2  n+l)Rf(tn,n) 


^„,  (2  m+l)iTZ     ^„,   (2  n+1)Try       ^,       . 
cos  -^ — 2h    —  "^ vj  ~  9vni»n) 


(IV-17) 
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IV-A.2   Concentration  Distribution 

The  form  of  the  transport  diffusion  equation  used  to  solve  for 
c"(y,z,t)  is  identical  to  (11-27),  or 

3c:-  _  -  Sfcl.  K  !!£:=  -  u"|i-  (IV-18) 

where  u"(y,z,t)  is  given  by  (IV-7).  It  is  assumed  that  c"(y,z,t)  is  of 
the  form 

c"(y,z.t)  =  c^  (y.z)  e''^*  (IV-19) 

where  c"(y,z)  is  complex  following  the  same  approach  used  with  u"(y,z,t). 
Substituting  (IV-19)  into  (IV-18)  and  requiring  zero  diffusive  flux  across 
the  boundaries  of  the  flow,  the  boundary  value  problem  is  then  stated  as 

8c" 
BC's:  a/  I  =  0 
^^   z  =  ±h 


3c" 


W 


I  =0 
y  =  ±w/2 


Mow.  c"  must  also  be  an  even  function  in  y  and  z  based  on  physical 
s 

arguments  and  is  therefore  assumed  to  be  of  the  form 

1=0    k=0 
Equation  (IV-21)  is  seen  to  satisfy  the  conditions  of  zero  diffusive  flux 
across  the  flow  boundaries.  Substituting  for  c^  using  (IV-21)  and  for 
u"  using  (IV-4)  and  (IV-7),  (IV-20)  becomes 
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00  00 


Yl  X]  '"^z  (f  )^  '  S  '^'^*  ^"'  ^ik  "=  T  "^  ^  = 


5,  =  0    k  =  0 


00  CO 


-  (|  Z  Z]  '™n  'COS  ii^iili^  cos  i^!^    -  g(™.n,> 
m  =  0  n  =  0  (IV-22) 


where  a      is  given  by  (IV-9).     The  unknown  Fourier  coefficients  a'^^  are 
solved  for  in  the  same  manner  as  was  used  to  determine  a^^.     The  detailed 
mathematics  of  the  solution  are  presented  in  Appendix  B,  the  result  of 
which  is  given  by 


(p+q) 


■M 


[(2  p.l)^-(2£)^][(2  q.l)^-(2k)^]re^[(^|i:^]^+e/^^f^]^  ^  ia} 

Thus,  the  solution  for  c"(y,z,t)  includes  four  infinite  summations 
arising  from  the  non-orthogonal  nature  of  the  cosine  functions 
representing  u^(y,z)  and  c^(y,z).  Again  introducing  Reynolds  analogy 
and  the  expression  for  the  vertical  and  transverse  mixing  times,  the 
solution  for  c"(y,z,t)  in  polar  form  as  a  function  of  the  pressure 
gradient  is 


CX)  OD  DO  UU 


£  =  0k  =  0p  =  0q  =  0 


/   ix-fJi+k)         l-nz  2k7Ty       iO, 

V^-^)  ^cos -y- cos  ^  e  ^  (IV-24) 

[(2  p+l)^-(2il)2][(2  q+l)2-(2k)2]R^^(p,q)R2^(£,k) 
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where, 

8^  =  at  -  ai(p,q)  -  az(Jl,k)  (IV-25) 

Ri(p,q)  is  given  by  (IV-12)  with  p  replacing  m  and  q 
replacing  n. 

^  'cz         cy 
ai(p,q)  is  given  in  (IV-13)  with  p  replacing  m  and  q 
replacing  n. 

a.(^k)  -   tan-'  (  ^   ,  ,.^    I   „  J  '"""' 

0,  k  =  iJ-  =  0 
^ilk  =  I  is,  k  =  0  or  £  =  0  (IV-28) 

1,  k  ?<  0,£?«  0 

The  factorv^.  is  necessarily  included  in  (IV-24)  to  properly  account 
for  the  a'  and  a^Q  terms,  and  to  insure  that  c"(y,z,t)  has  a  zero  mean 
over  the  flow  cross-sectional  area.  Detailed  mathematics  for  arriving 
at  (IV-24)  are  given  in  Appendix  B. 

The  expression  of  c"(y,z,t)  as  a  function  of  u^^^  is  determined  by 
substitution  for  K  in  (IV-24)  using  (IV-15)  to  yield. 


00  00  00  00 


11/         a.\       /     9C\       max 
c"(y.z.t)  =  (-  ^)  ^^r-f 


1=0       k=0       p=0       q=0 


1.      (  ir^^'^'^cos  ^^cos  ^iSIIe'^^ 

^£k^"^^  ^°^    h    ^°^    w  (IV-29) 

[(2  p+l)2-(2£)^I(2  q+l)2-(2k)2]R^^(p,q)Rj(£,k) 
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IV-A.3   Dispersive  Mass  Transport 

The  time  mean  dispersive  mass  transport  averaged  over  one  period  of 
oscillation  for  the  full  rectangular  section  as  shown  in  Figure  11  is 

given  by 

.w/2  ^h  V 

R,Cu"(y.z.t)]-R  Cc"(y,z,t)]  dz  dy  V   (IV-30] 


The  integrations  in  (IV-30)  are  carried  out  in  Appendix  B  and  it  is  noted 
that  the  non-orthogonality  of  the  cosine  functions  representing  u"  and  c", 
arising  from  the  form  of  their  respective  arguments,  once  again  produces 
two  additional  infinite  summations.  Thus,  using  (IV-11)  and  (IV-24)  for 
u"  and  c"  the  following  expression  for  the  longitudinal  dispersive  mass 
transport  over  the  region  -h  <  z  <  h  and  -w/2  <  y  <  w/2  is  determined 


CO  00 


^>,M-t)?^^^^^'Z  £  Z  Z  Z  Z 

i  =  0     k  =  0  m  =  0  n  =  0  p  =  0  q  =  0 


6i(Ji,k,m,n,p,q) 


where. 


X  =  ai(m.n)  -  ai(p,q)  (IV-32) 


6i(ii,k,m,n,p,q)  =  [(2p+l)^-(2£)^][(2q+l)^-(2k)^][(2n+l)^-  (2k)^]  • 

[(2m+l)2-(2£)2]  R'^^(m,n)R'f(p,q)R2(£,k)        (^^-33) 
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Now  the  dispersive  mass  transport  associated  with  the  open  rectangu- 
lar channel  comprising  the  lower  half  of  the  full  box  is  i5<m>j.  Therefore, 
the  Fickian  flux  representation  of  the  dispersive  mass  transport  in  the 
open  rectangular  channel  is  simply 

-<E^>^  (wh)  II  =  h<m>^  (IV-34) 

Replacing  <m>^  in  (IV-34)  with  the  right  hand  side  of  (IV-31)  and  solving 
for  <E^>-[-  yields 

00  00  00  00  OO  00 

a  =  0     k  =  0     m  =  0     n  =  0     p  =  0     q  =  0 


y     {(R2(£.k)-1)^^  cos  X  +  sin  X} 
^"^ _ — . (Iv-35) 

6i(£,  k,m,n,p,q) 

The  longitudinal  dispersion  coefficient  is  expressed  as  a  function 

of  u   by  substituting  for  K  in  (IV-35)  using  (IV-15)  to  give 
max 

00  00  OO  OO  00  OO 

^  '^         X.  =  0k  =  0m  =  0n  =  0p  =  0q  =  0 


yoJ(R2(Ji,k)-l)^  cos  A  +  sin  A} 

^k         -_ (IV-36) 

6i(£,k,m,n,p,q) 

To  facilitate  the  presentation  of  results  and  the  application  of 
(IV-36)  for  predictive  purposes  a  non-dimensional  relative  mixing  time. 


T',  is  introduced  as 
c 


1'==^=^  (IV-37) 

"  \  ^cy 
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where. 


T.  =  Zc^ 
'y   T 


(IV-38) 


Thus,  as  shown  by  (IV-37),  T'  is  a  measure  of  the  relative  effects  of 

vertical  and  lateral  shear  for  a  given  width-to-depth  ratio,  and  vertical 

and  lateral  eddy  diffusivities.  With  the  aid  of  (IV-37)  the  following 
definitions  can  then  be  made: 

r,(m.n)  =  {2uT')^R,(m,n)  =  {[^^1^]%  T^[i2n±DjL] 2^2,(2^1;)^    (IV-39) 


r,{Z,k)  =   (2TTT')2R,(£,k)  =  {(£TT)2+T'(kTT)2}2+(2TrT;)2 


;iV-40) 


ai(m,n)  =  tan   , 


2ttT: 


^'  = 


W^ 


j-(2m+1)TT-|2,-p,|-(2n+1)^-|2 

(-l)"^^"[ri(m.n)-(2uT;)^t^ 
_,  (2m+l)(2n+l)rj(m,n) 


m  =  0  n  =  0 


(IV-41) 


00     00 


+  <   2-nT. 


ill 


m+n 


J  (2m+lj(2n+l)ri(m,n) 


m  =  0  n  =  0 


(IV-42) 


These  relationships  are  then  applied  to  (IV-36)  to  obtain  an  expression 
for  the  non-dimensional  longitudinal  dispersion  coefficient,  E^,  defined 
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by  (111-58).  This  expression  represents  the  final  form  of  the  solution 
for  the  case  of  three  dimensional  oscillatory  shear  flow  and  is  there- 
fore presented  along  with  a  summary  of  the  associated  terms  as  follows: 


<E  >T      32T; 
p,   _       X  T  _  z 

^^^  jl=0k  =  0m=0n  =  0p  =  0q  =  0 


oo 

oo 

oo 

<xj                     oo                     00 

ZZZEZX 

y^k^[r2(£,k)-(2TTT;)2f^  cos   X  +  (2uT;)  sin  X} 


{IV-43) 


62(ii.k,m,n,p,q) 
where, 
62{il,k,m,n,p,q)   =  (ZirV  )'*6i(£,k,m,n,p,q) 

=  [(2p+l)^-(2il)2][(2q+l)2-(2k)2][(2m+l)2-(2£)^]   • 

[(2n+l)2-(2k)2]r^^(m,n)r^^(p,q)r2(Jl,k)  (IV-44) 

ri(m,n)  ={[-^^^1^+  r^li^L^y}^  +  (2uT;)2 
ri(p,q)  =  ri(m,n)  with  p  replacing  m  and  q  replacing  n 

r2(Ji.k)  ={(£TT)^+T;.(kw)^}^  +  (2ttT;)^ 

2TrT' 
ai(m,n)  =  tan'^j ^(^^^^ ^^^^^^^ ^(^^^^ )^^^  } 


ai(p,q)  =  ai(m,n)  with  p  replacing  m  and  q  replacing  n 
-   (-l)''^"[ri(m,n)-(2TrT;)2f^      ] 

2  1        V  \  ^ 


^     M   ^     ■    y     ■         (2m+l)(2n+l)ri(m,n)  f 

=  0     n  =  0  J 

(00  OQ 

^''"^z    ^  „4-7^     (2m+1)(2n+l)r,(m,n 
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'ik 


X  =  ai(ni,n)  -  ai(p,q) 

1  ,  k  ?f  0.  ^  ?*  0 
h  ,  k=Oor£=0 
0  ,  k  =  a  =  0 


The  syrmetry  of  this  solution  can  now  be  seen.  If,  for  example, 

T'  had  been  inversely  defined  as  l^yH^z'   *^^"  ^^^"^^^  ''°''^''  ^^''^ 
reduced  to  (IV-43)  with  T'  replacing  T^  everywhere,  and  the  newly  defined 
T'  multiplying  all  complementary  terms  in  ri,  rz,  and  aj.  This  form  of 
{IV-43)  would  predict  the  same  value  for  E^  as  the  original  formulation 

provided  that  u   and  T  were  the  same  and  that  the  new  T'  equaled  the 
"  max  ■^ 

old  T'  and  the  new  T'  equaled  the  old  T'   Stated  another  way,  for  the 
z  ~^^~  c 

same  u   and  T,  a  channel  whose  half-width  was  twice  its  depth  would 
max 

produce  the  same  dispersive  mass  transport  as  a  channel  whose  depth  was 
twice  its  half-width  provided  that  the  ratio  yR^  for  the  first  case 
equaled  K  /R  in  the  second  case. 

The  functional  behavior  of  E^  as  given  by  (IV-43)  for  three  dimen- 
sional oscillatory  shear  flow  is  illustrated  by  Figure  12.  In  this 
figure  E'  is  plotted  against  the  non-dimensional  vertical  mixing  time, 
T',  for  fixed  values  of  the  relative  mixing  time,  TJ..  Because  the 
solution  for  E'  contains  six  infinite  summations  nested  in  series  some 
limitations  were  necessary  in  carrying  out  the  required  computations. 
By  varying  the  upper  bound  on  each  of  the  sums  it  was  demonstrated  that 
the  solution  converges  upward  to  its  limiting  value.  The  convergence 
occurs  reasonably  rapidly;  however,  above  an  upper  limit  of  10  terms  for 
each  sum  the  rate  of  convergence  is  slowed  considerably.  The  results 
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FIGURE  12     NON-DIMENSIONAL  LONGITUDINAL  DISPERSION  COEFFICIENT 

E',  AS  A  FUNCTION  OF  T'   FOR  FULL  RANGE  OF  T'   VALUES 
X  z  c 
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presented  in  Figure  12  and  in  subsequent  figures  were  computed  using  an 
upper  bound  of  10  on  each  sum;  therefore,  each  value  of  E^  includes 
approximately  1.2  x  10^  terms.  An  estimate  of  the  error  between  the  com- 
puted values  of  E'  and  the  values  predicted  by  the  complete  series 
represented  by  (IV-43)  is  obtained  by  comparing  the  limiting  curve  for  a 
very  wide  channel  in  Figure  12  with  the  corresponding  curve  for  an 
infinitely  wide  channel  shown  in  Figure  7.  The  limiting  curve  in  Figure 
12  is  identified  by  T'  <  1  x  10"^  on  the  left  hand  side  of  the  peak  and 
by  T'  <  0.001  on  the  right  hand  side  of  the  peak.  A  comparison  of  the  two 
curves  suggests  that  for  T'  <  1  x  10"  the  solution  for  E'  given  by 
{IV-43)  ^ery   nearly  approximates  for  all  T^  the  solution  given  by  (III-59) 
for  E'  in  an  infintely  wide  channel.  The  general  shapes  of  the  two  curves 
are  nearly  identical  with  both  peaks  occurring  at  T^  =  1.58.  For  the 
infinitely  wide  channel  the  peak  value  of  E^  is  3.27  x  10'^  whereas  for 
the  limiting  curve  in  Figure  12  it  is  3.07  x  10"^  or  6  per  cent  below 
the  infinitely  wide  case.  Thus,  it  appears  that  the  three  dimensional 
solution  for  E'  approaches  the  two  dimensional  solution  for  E'  in  the 

X  X 

limit  as  T'  -^  0,  and  it  is  estimated  that  values  predicted  for  E'  in 
c  X 

Figure  12  are  low  but  generally  within  10  per  cent  of  their  actual 
values  as  determined  by  (I\/-43).  This  does  not  apply  to  the  range 
T'  >  100  where  the  Figure  12  values  are  as  much  as  70  per  cent  low. 
However,  this  is  not  a  region  of  practical  significance  since  values 
of  T'  in  nature  are  generally  <100. 

To  aid  in  the  interpretation  of  results,  the  eight  curves  shown  in 
Figure  12  have  been  separated  into  two  groups  and  are  replotted  in 
Figures  13  and  14.  Figure  13  shows  the  E'  resonance  curves  for  T^  =  1.0, 
0.67,  0.33,  and  0.1  while  Figure  14  shows  the  corresponding  curves  for 
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FIGURE  13  NON-DIMENSIONAL  LONGITUDINAL  DISPERSION  COEFFICIENT 
E'  AS  A  FUNCTION  OF  T^  FOR  T^  =  1.0,  0.67,  0.33 

AND  0.1 
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FIGURE  14     NON-DIMENSIONAL  LONGITUDINAL  DISPERSION  COEFFICIENT 
E'  AS  A  FUNCTION  OF  T^  FOR  T"   =  0.1,  0.01,  0.001, 

0.0001,  and  <1  x  10"6 
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T'  =  0.1,  0.01,  0.001,  0.0001,  and  <  1  X  10"  .  Note  that  the  curve  for 

T'  =0.1  is  repeated  to  provide  a  common  reference.  If  for  the  purposes 

of  interpretation  it  is  assumed  that  K  =  K  then  T'  may  be  written  as 

z    y      c 


"^c  =  ^W^'^' 


(IV-45) 


which  is  simply  the  inverse  square  of  the  channel  half-width  to  depth  ratio. 

Referring  to  Figures  12,  13,  and  14  with  the  aid  of  (IV-45)  it  is  seen 

that  as  the  half-width  of  the  channel  is  increased  from  a  value  equal  to 

the  depth  the  resonant  peak  shifts  from  T'  =  4.5  leftward  toward  T^  =  1.58 

at  which  it  occurs  when  T'  <  0.001.  Concurrent  with  the  shift  is  an  ini- 

c  - 

tial  increase  in  the  peak  value  of  E'  as  T'  is  decreased  from  1.0  to  0.25 
followed  by  a  decrease  in  the  peak  value  as  T'  is  further  decreased  to 
0.001.  Decreasing  T'  below  this  value  has  no  effect  on  E^  for  T^  >  1.2. 
A  summary  of  the  maximum  value  of  E'  achieved  and  the  value  of  V  at  which 
it  occurs  for  the  corresponding  T'  is  given  in  Table  1.  These  data  are 
also  presented  graphically  in  Figures  15  and  16. 


TABLE  1 

PREDICTED  RESONANT  PEAK 

DATA  FOR 

LONGITUDINAL  DISPERSION 

COEFFICIENT 

IN 

THREE  DIMENSIONAL  OSCILLATORY  SHEAR  FLOW 

Non-Dimensional 

Peak  Value 

Vertical  Mixing 

of  Non-Dimensional 

Relative  Mixing 

Time,  T'  at  which 

Longitudinal  Dispersion 

Time,  T" 

z 

Coefficient,  E'  x  10^ 

Maximum  Occurs 

1.00 

4.50 

339 

0.67 

3.75 

340 

0.40 

2.40 

350 

0.33 

2.20 

353 

0.25 

1.95 

354 

0.20 

1.80 

352.5 

0.10 

1.65 

346 

0.01 

1.60 

325 

0.001 

1.58 

310 

0.0001 

1.58 

307 

0.00001 

1.58 

307 

0.000001 

1.58   ~ 

307 
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FIGURE  15     NON-DIMENSIONAL  VERTICAL  MIXING  TIME,  T^,  AT 

WHICH  MAXIMUM  E'   OCCURS  AS  A  FUNCTION  OF  T' 
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FIGURE   16     PEAK  E^  AS  A  FUNCTION  OF  T^ 
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The  effect  of  T'  on  the  general  shape  of  the  resonance  curve  for 
E'  is  complex.  To  understand  this  behavior  it  is  best  to  treat  the 
right  and  left  hand  sides  of  each  curve  separately  and  to  note  that  for 
the  two  extreme  channel  geometries  corresponding  to  T^  =  1.0  and  T^  < 
1  X  10"^  the  left  hand  sides  of  the  respective  curves  are  coincident 
up  to  T'  =  1.0.  Looking  at  the  right  hand  side  of  the  curves  in 
Figure  12  it  is  seen  that  as  T'  is  decreased  from  1.0  to  0.001  the  curves 
are  depressed  toward  the  limiting  two  dimensional  curve,  merging  with 
it  at  correspondingly  lower  values  of  T'  as  T'  is  decreased. 

For  values  of  T'  <  0.001  the  right  hand  side  of  the  curve  remains 

unchanged.  Thus,  for  values  of  T^  >  7  the  effect  of  lateral  shear  on 

longitudinal  dispersive  mass  transport  in  oscillatory  flow  is  greatest 

for  width-to-depth  ratios  corresponding  to  T^  =  1.0  and  decreases  as 

the  width-to-depth  ratio  is  increased  to  a  value  corresponding  to  T^  = 

0.001.  For  T'  <  0.001  the  lateral  shear  has  no  effect  and  the  dispersive 
c 

mechanism  behaves  as  a  two  dimensional  shear  flow  phenomenon. 

The  behavior  of  the  left  hand  side  of  the  E^  resonance  curves  is 
most  easily  seen  from  Figures  13  and  14.  As  T^  is  initially  decreased 
from  1.0  to  0.33  the  curve  shifts  to  the  left  thereby  predicting  higher 
values  of  E'  for  corresponding  values  of  T'   Decreasing  T'  beyond  this 

X  Z  L 

point,  however,  causes  the  upper  portion  of  the  curve  to  shift  back 
toward  the  right  predicting  lower  values  of  E'  while  the  lower  portion 
of  the  curve  is  flattened  predicting  higher  values  of  E'  for  small 
values  of  T'.  Referring  to  Figure  14,  as  T'  is  decreased  from  0.1  to 
1  X  10"^  the  upper  portion  of  the  curve  continues  to  be  supressed 
toward  the  right  and  merges  with  the  limiting  two  dimensional  curve  at 
correspondingly  lower  values  of  T'  while  the  lower  extreme  portion  of 
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the  curve  continues  to  flatten  until  it  also  merges  with  the  limiting 
two  dimensional   curve.     For  TJ.  <  1  x  10"     the  left  hand  side  remains 
unchanged. 

It  is  observed  that  provided  one  knows  the  values  of  the  non- 
dimensional  vertical  and  lateral  mixing  times,  T^  and  T^,  it  is  possible 
to  predict  the  limiting  width-to-depth  ratio  for  which  the  three  dimen- 
sional oscillatory  flow  case  reduces  to  the  more  easily  applied  tvo 
dimensional  case.     The  general  condition  for  which  this  is  true  may  be 
stated  as 

T'  <  T'  (IV-46) 

c  -     o» 

where  T'^  is  the  limiting  value  of  T^  for  which  the  two  dimensional 

oscillatory  dispersion  solution  is  valid.  Using  the  definitions  of 

T'.  T  ,  and  T   (IV-46)  can  then  be  written  as 
c'  cz'     cy  ' 

w   K 
z 
or, 

{IV-47) 


It  is  noted  from  (IV-47)  that  increases  in  K  due  to  secondary 
currents  and  horizontal  eddies  produced  by  channel  irregularities  will 
increase  the  value  of  w/h  for  which  the  two  dimensional  solution  applies 
Similarly,  decreases  in  K  due  to  density  stratification  will  have  the 
same  effect.  Values  of  T'^  are  obtained  from  Figure  12.  It  is  assumed 
that  the  channel  is  straight  and  that  the  fluid  is  homogeneous  with  no 
density  stratification.  For  these  conditions  Ward  (23)  has  shown  in  a 
detailed  set  of  laboratory  experiments  that  the  time  averaged  ratio 
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1/ 
of  -X  in  an  oscillatory  flow  is  approximately  twice  the  value  determined 

by  Elder  (4)  for  unidirectional  flow.     Using  the  mean  of  four  experimental 
runs  from  Ward's  data  it  is  assumed  that 

JL=6.6  (IV-48) 

The  predicted  limiting  width-to-depth  ratios  based  upon  these  assumptions 
are  presented  in  Table  2. 

TABLE  2  LIMITING  WIDTH-TO-DEPTH  RATIOS  FOR 
APPLICATION  OF  THE  OSCILLATORY  TWO 
DIMENSIONAL  SHEAR  FLOW  DISPERSION 
COEFFICIENT 


T'  values 
which  val 

for 
id 

T'  values  for 
c 

which  valid 

Width-to-depth 
Ratio,  (w/h) 

>6.0 

<0.01 

>51 

>1.2 

<0.001 

>162 

>0.35 

<G.0001 

>514 

>0.04 

<0. 00001 

>1625 

+00 

<0. 000001 

>5138 

<1.2 

1.0 

5.1 

For  given  values  of  T'  and  T'  Table  2  is  used  to  determine  if  (w/h)  is 
^  z     c 

large  enough  to  allow  the  use  of  the  more  simple  solution  for  E^  pre- 
sented in  III-B. 

Before  passing  from  the  present  section  it  is  re-emphasized  that 
while  the  preceding  remarks  have  been  directed  toward  channel  geometries 
in  which  the  half-width  is  always  equal  to  or  greater  than  the  depth, 
they  apply  equally  to  channel  geometries  in  which  the  depth  is  always 
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equal  to  or  greater  than  the  half -width  provided  the  conditions  for 

syirmetry  previously  discussed  are  met.  Thus,  Figures  12,  13,  and  14  could 

be  used  equally  well  for  cases  where  |^  <1  provided  T^  were  replaced 

with  T'  and  T'  were  redefined  as  T'/T'.  Perhaps  the  best  way  to  show 
y     c  y  z 

the  syimetry  of  the  solution  about  T'  =  1  is  by  means  of  a  plot  such  as 
Figure  17.  In  this  figure  isolines  of  E^  for  channel  geometries  skewed 
both  in  depth  and  width  are  plotted.  The  lower  half  of  this  figure 
corresponds  to  geometries  in  which  ^  >  1  whereas  the  upper  half 
corresponds  to  geometries  in  which  2h"-  ^'  provided  K^  =  K  . 
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9    Maximum     Value     of     E^  =  354  X  10"^  occuring    at    T(,'=0.25 
and     T^  =  1.95 


FIGURE   17     SYMMETRIC  DISPLAY  OF   ISOLINES  OF  E^  x  10 

IN  THREE  DIMENSIONAL  OSCILLATORY  SHEAR  FLOW 
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IV-B  Steady  Unidirectional  Flow 
The  following  analysis  is  simplified  by  using,  where  applicable, 
the  results  of  the  previous  section. 
IV-B.l   Velocity  Distribution 

The  coordinate  system  used  and  the  boundary  value  problem  formula- 
tion for  the  steady  flow  case  are  identical  to  that  described  in  IV-A 
except  for  the  exclusion  of  the  local  acceleration  term.  Thus,  for  no 
slip  conditions  on  the  perimeter  of  the  full  rectangular  section  and  a 
constant  pressure  gradient  as  defined  by  (III-2)  the  flow  boundary  value 
problem  for  u(y,z)  is 

PDF.    -  8iii+  e  9!ii  =  K  (IV-49) 

1-Hl-    ^z  8z2   ^y  3y2 

BC's:   u(±w/2,z)  =  0 
u(y,±h)   =  0 

Proceeding  in  the  same  manner  as  before,  assume  a  solution  of  the  form 

OO  00 

>           >                       (2  m+l)iTZ       ^   (2  n+l)TTy  ,,,/  rn\ 

u(y,z)  =  ^  ^  a^^  cos        g^         "'    ^  ^^^"^^^ 

m  =  0  n  =  0 

which  satisfies  the  boundary  conditions  and  provides  for  zero  vertical 
shear  stress  at  the  free  surface  and  zero  lateral  shear  stress  at  the 
centerline  of  the  channel.  Equation  (IV-50)  is  then  substituted  into 
(IV-49)  to  yield 
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00  °° 


m'  =  0    n  =  0 


,      cos  <i^i2^    cos  ilJili2L    ,  .K 
mn  Zh  w 


(IV-51) 


A  comparison  of  (IV-51)  with  (B-3)  in  Appendix  B  shows  them  to  be 
of  the  same  general  form.  Therefore,  the  solutions  for  u(y,z)  and  u"(y,z) 
may  be  written  directly  as 


00  00 


-16K(-ir"coslg!!|l^cosl^aalSl    ^ 


and. 


00  oo 


■■(^•^'  =  2^  Yl  ^™n  f  ™=  ^^S^    -=  ^^'^  -g(m,n))  (IV-53) 


m  =  0  n  =  0 
where , 


3       _    -16K(-ir-^"  (IV.54) 

">"        .^(2m+l)(2n+l){e^[i^ni±lJZ]a,  E^[i^I!±II^} 

The  solution  for  u"(y,z)  as  a  function  of  the  pressure  gradient  is 
put  in  its  final  form  by  utilizing  the  Reynolds  analogy  and  introducing 
the  mixing  times  T^^,  T  ,  and  T^  where  in  this  problem  T^  can  only  be 
written  as 
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T'  =  =^  (IV-55) 

since  there  is  no  period  of  oscillation.  The  resulting  expression  is 


^    m=0  n=0  """  ^ 

(2m+l)TTZ     (2n+l)Try  „,  „>  ( 
cos  ^  2h    ^°'    w   -Q^'"'")  j        (IV-56) 

where, 

R   =  {C(2"i^l)^]%  T'  [i^^]'}       (IV-57) 
mn      ^      c    t 

To  express  u"  as  a  function  of  u^^^, (IV-52)  is  evaluated  at 
y  =  z  =  0  and  the  result  solved  for  K  to  yield 

2 

K=.^^f>L^  (IV-58) 

where. 


CO        00 

^m+n 


m=0  n=0 
Substituting  (IV-58)  into  (IV-56)  yields 


00  CO 

iHi+n 


..,       ,       "max    V  V"  .     (-ir "   . 
u"(y.z)  =  ^T-    ^   ^  (2m+l)(2n+l)R^^ 

m=0     n=0  L 

cos  ^2^^    COS  -^-^^^  -  g(m.n)    j  (IV-60) 
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IV-B.2   Concentration  Distribution 

The  boundary  value  problem  for  the  steady  state  solution  of  c"(y,z) 
using  a  no  flux  condition  on  the  perimeter  of  the  section  is 


PDE: 


K 


9^c' 


+  K. 


.9^c"  _  ..„  9c 


Z  9z2   'y  9y2 


=  u 


9C 


(IV-61) 


BC's: 


9c" 

9Z 


=  0 
z  =  ±h 


9c^ 

9y 


=  0 
y  =  ±w/2 


Then  assuming  a  solution  to  (IV-61)  of  the  form 
c"(y.z)  = 


00  00 


1  =  0     k  =  0 
and  substituting  into  (IV-61),  the  result  is 


z I  z I   z  n      y      VI 


*£k 


S,  =  0  k  =  0 


(--) 
^  9^^ 


00  00 


m  =  0  n  =  0 


IV-62) 


.  «mn["^  -^^^  ^°^  ^^^  -9(-'")>   (^V-63) 


Equation  (IV-63)  is  of  the  same  form  as  (IV-22)  thereby  allowing 
the  solution  for  c"(y,z)  to  be  written  as 

^  'u,,(-l)-<^^^)cos^cos2M 


OO  OO  CO 


c"(y.z)  =  (|f)  ^ 


'Sik' 


63(^1, k,p,q] 


1  =  0     k  =  0p  =  0q  =  0 


(IV-64) 
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where, 

«^3(Jl.k.p.q)  =  [(2p+l)^-(2£)^][(2q+l)2-(2k)2]  • 

a^li^2^y^  l^lila^yn'K^i^V^  l^i^yy  (IV.65) 

1    k  ?«  0,  Ji  ?«  0 
^nj^   \  h  k  =  0  or  ji  =  0 

0    k  =  Jl  =  0 

To  express  c"  as  a  function  of  T   and  T'  Reynolds  analogy  is  applied 
as  before  and  the  mixing  time  scales  are  introduced  so  that  (IV-64)  becomes 
..  256KT^    «>     CO     00     CO 


i  =  0     k  =  0  p  =  0  q  =  0 


.,,(-1)-"-^)  cos  ^  cos- 


2kTr\ 


(IV-66) 


6i»{i,k,p,q) 
where, 

<S.(5<,k,p,q)  =  [(2p+l)2-(2£)2][(2q+l)2-(2k)2]  R^^  Rp^        (IV-67) 

hk   "  ^^^^^'  ■"  "^c  ^'"'^^'^  ^^^-^^^ 


The  concentration  distribution  is  expressed  as  a  function  of  u   by 

^  max  "^ 


applying  (IV-58)  to  (IV-66)  to  yield 
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00  CO  00  00 

^ '    ^ '     \  I    ^ 1 


9C'       ^2 


TT^    ({>' 


l=Q       k=0       p=0       q=C 


^Zk{-1)'^^^^'    cos-[p  cos    w  (IV-70) 

6j5,,k,p,q) 

IV-B.3   Dispersive  Mass  Transport 

The  dispersive  mass  transport  resulting  from  a  unidirectional  flow 
through  the  full  rectangular  section  of  Figure  11  is  given  by 

m  =   \     \    u"(y,z)  c"(y,z)  dy  dz      (IV-71) 
-'-h    -w/2 

Using  (IV-56)  for  u"(y,z)  and  (IV-66)  for  c"(y,z)  the  integrations 
in  (IV-71)  are  carried  out  in  the  same  manner  as  described  in  Appendix 
B.  Then,  as  discussed  in  IV-A  the  open  channel  dispersive  mass  trans- 
port equals  one-half  m  and  is  described  as  a  Fickian  flux  by 

-E^  (wh)  ^=  hfn  (IV-72) 

The  resulting  expression  for  E.  as  a  function  of  the  pressure  gradient  is 

00  00  00  00  CO  oo 

1=0     k=0    m=0     n=0     p=n     q=0 
where, 

63{£,k,m,n,p,q)  =  [(2p+l)2-(2il)2l(2q+l)2-(2k)2:  • 
C(2ni+l)2-(2£)^X(2n+l)2-(2k)2:  R^^  Rpq  R^,^  (IV-74) 
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The  longitudinal  dispersion  coefficient  is  expressed  as  a  function  of 

u   by  substituting  for  K  in  (IV-73)  using  (IV-58).  This  expression 
max  -^ 

is  put  in  non-dimensional  form  by  defining 


"^L   u2   T 
max  cz 


{IV-75) 


The  final  form  of  the  solution  is  summarized  below. 
E, 


OO  00  CO  oo  oo 


64 


/■  y.  ZZ  Z^  Z^    Z^  65(il,k,m.n!p.q)    ^^^-^S) 


n.      _  -L 

L         "may   ^r7        ^^^V^ 

where, 

65   (il,k,m.n,p,q)   =  [(2p+l  )2-(2£)2][(2q+l  )2-(2k)2][(2m+l  )2-(2£)n- 
[(2n.l)^-(2k)^]R^^RpqR,, 


R      =  R      with  p  replacing  m  and  q  replacing  n 
pq        mn  ^       r  » 


hk  =  t^(^^)'  ■"  "^c  ^""^^'^ 


*'     =  Z_.   y     ■  T2iiTT7T 


m+n 


2n+l)R 


mn 


'«,k 


m  =  0     n  =  0 

'"I  k  ?«  0,  £  ?«  0 

%  k  =  0,  or  £  =  0 

0  k  =  £  =  0 
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It  should  be  noted  that  the  above  solution  is  symmetrical  about  T'  =  1 
in  exactly  the  same  manner  as  the  three  dimensional  oscillatory  shear 
flow  solution,  and  although  the  discussion  to  follow  is  restricted  to 
channels  whose  half-widths  are  equal  to  or  greater  than  their  depths  the 
same  remarks  apply  to  those  channels  whose  geometries  are  skewed  in  the 
opposite  di recti oa. 

Equation  (IV-76)  is  plotted  in  Figure  18  vs.  T'.  In  performing  the 
computations  a  summation  limit  of  12  was  used  on  each  of  the  six  infinite 
summations  given  by  (IV-76).  By  adjusting  this  value  in  a  manner  similar 
to  that  discussed  for  the  oscillatory  case  it  was  concluded  that  the 
solution  as  presented  in  Figure  18  is  ^ery   near  that  given  by  the  infinite 
summations. 

It  is  seen  from  Figure  18  that  for  a  unidirectional  shear  flow  in  a 
rectangular  channel  of  finite  width  the  longitudinal  dispersion  coeffi- 
cient for  wide  channels  becomes  very  large.  By  comparison,  if  the 
definition  of  Ej"  given  by  (IV-75)  is  applied  to  the  two  dimensional 
solution  for  E.  given  by  (III-18)  the  result  is 

E|^  =  gig  =  constant  (IV-77) 

which  is  shown  in  Figure  18  by  the  dotted  line.  The  difference  between 
the  behavior  of  the  two  solutions  is  dramatic  and  for  unidirectional 
flows  supports  Fischer's  theory  (6)  regarding  the  dominant  effect  of 
lateral  shear  on  longitudinal  dispersion  in  wide  channels.  The  fact  that 
the  three  dimensional  solution  does  not  approach  the  two  dimensional 
solution  in  the  limit  as  T'  ^  0  is  not  surprising  if  one  realizes  that 
the  two  problems  are  completely  separate  and  independent  of  one  another. 
The  three  dimensional  problem  has  by  definition  a  lateral  shear  effect 


98 

10^ 

\ 

10 

\ 

\ 

-1 

UJ 

1 

K 
O 

J^  lo" 

\.^^^^^ 

10^ 

—  ■ 

E;J_2  DIMENSIONAL)  =  8.46X10"^         ^^^^^ 

-3 

,        1        1        1        1  1 

10 

1 

0^ 

i6'    lo''    10"^    10^    lo'     1 

t;  ■ 

FIGURE  18  NON-DIMENSIONAL  LONGITUDINAL  DISPERSION  COEFFICIENT 
E^,  AS  A  FUNCTION  OF  T;.  FOR  THREE  DIMENSIONAL 

UNIDIRECTIONAL  SHEAR  FLOW 

99 


which  no  matter  how  wide  the  channel,  is  forever  present.   Conversely, 
the  two  dimensional  problem  by  definition  has  no  such  effect.  The 
presence  of  this  shear  coupled  with  a  unidirectional  flow  and  a  long 
time  scale  of  lateral  mixing  must  produce  a  significantly  higher 
longitudinal  dispersive  mass  transport  than  would  be  produced  in  the 
case  where  the  lateral  effect  is  non-existent. 


CHAPTER  V 
COMPARISON  OF  RESULTS  WITH  FIELD  AND.. 
LABORATORY  DATA 

V-A  Unidirectional  Flow 

As  an  evaluation  of  the  predictive  reliability  of  the  solution  for  the 
longitudinal  dispersion  coefficient,  E,  ,  in  three  dimensional  unidirectional 
shear  flow,  (IV-76)  is  compared  with  field  and  laboratory  data  as  analyzed 
by  six  investigators:  Fischer  (6,  8,  13),  Miller  and  Richardson  (24), 
Godfrey  and  Frederick  (25),  and  Thackston  and  Krenkel  (26).  As  shown  in 
Figure  19  the  data  are  very  scattered,  and  based  on  first  impressions  it 
can  only  be  said  that  the  solution  correctly  predicts  an  increasing  E' 
with  a  decreasing  T".  However,  the  assumptions  made  in  plotting  the  data 
are  many  and  the  conditions  for  which  they  were  obtained,  particularly  in 
the  field,  vary  greatly  from  the  assumed  straight  rectangular  geometry 
used  in  the  analytical  solution.  Values  for  the  longitudinal  dispersion 
coefficient  cannot  be  measured  directly  and  are  obtained  only  after 
considerable  reduction  of  data  which  necessarily  involves  subjective 
judgments.  This  fact  is  illustrated  by  comparing  the  observed  E'  values  in 
Figure  19  as  determined  by  Fischer  (8),  Godfrey  and  Frederick  (25),  and 
Thackston  and  Krenkel  (26).  All  of  these  values  were  reduced  from  the 
same  base  data. 

Perhaps  the  most  limiting  qualification  of  all  these  data  is  that  none 
include  measured  values  of  both  the  vertical  and  lateral  eddy  diffusivities. 
Only  one  set  of  data,  that  of  Miller  and  Richardson  (24),  includes  measured 
values  of  either  K  or  K  .  For  their  experiment  they  chose  to  measure  the 
lateral  eddy  diffusivity,  Ky.  The  lack  of  measured  K  and  K^  values  is 
unfortunate  since  it  forces  one  to  assume  values  for  T  ,  T  ,  and  T'. 
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Values  of  u   ,  the  maximum  cross-sectional  velocity,  are  also 
max 

needed  for  the  comparison.  Three  of  the  six  data  sets,  Fischer  (6,  8), 
and  Miller  and  Richardson  (24),  provide  cross-sectional  mean  velocity 
data  for  each  experimental  run  and  a  representative  isovel  distribution 
for  one  run.  Using  this  information  it  was  determined  that  the  relation- 
ship 

".ax  =  3/2  "  (V-l) 

predicted  by  the  theory  used  in  the  present  work  was  in  reasonable  agree- 
ment with  the  measured  data.  Additional  assumptions  were  used  to  obtain 

u   values  for  the  Fischer  (13)  data  and  these  will  be  discussed  separately, 
max 

Thus,  noting  the  above  considerations,  the  following  assumptions 
were  made  in  plotting  the  data  shown  in  Figure  19: 

(1)  ic^  =  0.063  u*h 

(2)  K  =  0.23  u^h  (except  where  measured) 

(3)  W3/2" 

The  expressions  used  to  estimate  K  and  K  are  based  on  the  results  of 
Elder  (4).  A  brief  description  of  the  data  follows. 

The  closed  circles  in  Figure  19  represent  data  obtained  by 
Fischer  (6)  in  a  detailed  set  of  three  dimensional  shear  flow  laboratory 
experiments.  The  flow  cross  section  was  trapezoidal  with  a  smooth  bottom 
and  rough  sides  achieved  by  carefully  stacking  an  embankment  of  1  inch 
stone  along  the  sides  of  the  flume.  Flow  depths  ranged  from  2.1  to  4.7 
cm  while  half-widths  ranged  from  9.5  to  21.5  cm.  Measured  values  of  E|^ 
used  to  compute  E/  are  those  given  by  Fischer  for  Series  2800  through 
3200  and  3400  runs  using  the  diffusive  transport  method.  The  degree  of 
correlation  between  these  data  and  the  predicted  values  of  E^^   using 
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(IV-76)  is  surprisingly  good  considering  the  assumptions  made  and  the 
trapezoidal  cross  section  of  the  flume. 

The  open  circle  data  in  Figure  19  were  reported  by  Fischer  (8)  who 
analyzed  field  data  collected  by  the  U.  S.  Geological  Survey.  A  detailed 
description  of  this  set  of  experiments  conducted  in  several  natural 
streams  is  given  by  Godfrey  and  Frederick  (25)  whose  results  are  shown 
in  Figure  19  by  the  closed  squares.  Thackston  and  Krenkel  (26)  also  used 
these  data  to  calculate  observed  values  of  E.  which  are  shown  by  the  open 
squares  in  Figure  19.  These  three  sets  of  data  are  widely  scattered 
with  most  points  falling  above  the  curve  predicted  by  (IV-76).  The 
general  lack  of  agreement  between  predicted  and  observed  values  of  E|^ 
for  these  data  is  not  surprising  in  view  of  the  assumptions  made  to 
determine  the  lateral  and  vertical  mixing  times  and  the  non-uniform 
nature  of  the  natural  channels.  The  largest  stream  included  in  the  data 
is  the  Clinch  River,  Tennessee  having  a  discharge  of  300  cfs  and  a  width- 
to-depth  ratio  of  nearly  60. 

Fischer  (13)  also  reports  several  experimentally  determined  values 

of  E.  as  determined  by  Thomas  (27),  Fischer  (28),  and  Yotsukura  et  al .  (29) 

for  the  Chicago  Ship  Canal,  Green-Duwamish  River,  and  Missouri  River 

respectively.  These  values  are  plotted  in  Figure  19  with  closed  triangles. 

To  obtain  values  of  u    for  these  data  it  was  assumed  that  the  cross- 
ma  x 

sectional  mean  velocity  and  the  shear  velocity  are  related  by 


PW  G  (v-2) 


where  f^^   is  the  Darcy-Weisbach  friction  factor.  Then  by  using  a  value 

of  0.08  for  fpij,  and  the  previously  assumed  relationship  between  u  and 

u„,„.  the  result  is 
max 
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"max  =  15  "*  (V-3) 

The  largest  channel  included  in  these  data  is  the  Missouri  River  having  a 
width-to-depth  ratio  of  74. 

The  final  set  of  data  shown  in  Figure  19  is  from  a  series  of 
laboratory  experiments  conducted  by  Miller  and  Richardson  (24).  These  data 
are  shown  by  open  triangles  in  the  figure  and  all  fall  along  a  narrow 
band  of  large  T'  values.  The  experiments  were  conducted  in  a  2  ft  wide 
by  60  ft  long  recirculating  flume  with  a  width-to-depth  ratio  of  approxi- 
mately 4.7.  Turbulent  flow  conditions  were  achieved  by  placing  artificial 
roughness  elements  at  the  bed.  Roughness  element  heights  ranged  in  three 
sizes  from  0.43  cm  to  2.86  cm  compared  to  a  nearly  constant  flow  depth 
of  approximately  13  cm.  Cross-sectional  mean  velocities  ranged  from 
30.5  to  77.4  cm/sec  with  the  higher  velocity  runs  using  the  largest 
roughness  elements.  Experimentally  determined  values  of  K  were  reported 
by  the  authors  and  were  used  here  to  determine  the  lateral  mixing  time. 
The  scatter  of  the  Miller  and  Richardson  data  is  considerable.  However, 
the  data  for  Runs  8  and  9  suggest  that  the  conditions  for  which  a 
longitudinal  dispersion  coefficient  may  be  defined  were  exceeded.  Mean 
velocities  for  these  two  runs  were  53.3  and  77.4  cm/sec  respectively  with 
flow  depths  of  12.6  and  13.1  cm,  and  roughness  element  heights  of  2.86 
cm.  Under  these  conditions  the  flow  would  begin  to  lose  its  orderly  shear 
profile  characteristics  and  would  certainly  exhibit  longitudinal  non- 
uniformities  due  to  the  relative  size  of  the  roughness  elements  compared 
to  the  depth  of  flow.  The  authors  cite  this  factor  in  addition  to  the 
entrapment  of  tracer  behind  the  roughness  elements  and  its  subsequent 
release  to  the  main  stream  as  possible  reasons  for  the  abnormally  high 
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values  of  E, .  In  either  case  it  seems  reasonable  to  place  a  reduced 
emphasis  on  the  data  from  Runs  8  and  9  which  correspond  to  the  two  top 
open  triangles  in  Figure  19. 

None  of  the  experimental  data  shov/n  correlate  particularly  well 
with  the  analytical  solution  obtained  in  Chapter  IV  for  the  dispersion 
coefficient.  However,  the  data  are  widely  scattered  with  no  well  defined 
trend  except  for  a  generally  increasing  E.'  with  decreasing  T'.  The 
importance  of  measuring  K  and  K  values  concurrent  with  E,  data  is 
definitely  emphasized  by  this  behavior.  In  view  of  the  idealized  nature 
of  the  solution  and  the  assumptions  necessary  for  a  comparison  to  be  made 
between  it  and  the  experimental  data,  it  is  felt  that  the  results  are 
reasonable.  Certainly,  (IV-76)  provides  a  general ,  predictive  expression 
for  the  longitudinal  dispersion  coefficient  as  a  function  of  both  vertical 
and  lateral  shear  and  diffusive  effects. 
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V-B  Oscillatory  Flow 

The  comparison  of  the  predicted  two  and  three  dimensional  oscillatory 
shear  flow  dispersion  coefficients  with  field  and  laboratory  data  is  also 
difficult  due  to  a  general  lack  of  adequate  data.  Because  of  the  sharply 
resonant  behavior  of  longitudinal  dispersion  in  oscillatory  flow  it  is 
essential  that  measured  values  for  K^  and,  where  applicable,  K  be  obtained 
concurrently  with  the  experimental  determination  of  <Ej^>y  Fortunately,  this 
has  been  done  by  Bowden  (5)  for  the  Mersey  Estuary,  Liverpool  Bay,  and  a 
region  off  the  Cumberland  Coast,  and  also  by  Bowles  et  al .  (30)  for  a  por- 
tion of  the  English  Channel  off  the  Dorset  Coast. 

Three  of  these  field  locations  may  be  regarded  as  two  dimensional 
oscillatory  shear  flow  situations  while  the  fourth,  being  the  Mersey 
Narrows,  must  be  treated  as  a  three  dimensional  oscillatory  flow.  Ob- 
served data  for  the  three  two  dimensional  cases  are  summarized  in  Table  3. 


TABLE  3  OBSERVED  DISPERSION  DATA  AND  PREDICTED 

<E  >,  VALUES  USING  (III-59) 
X  T 


LOCATION 

Parameters 

Liverpool  Bay 

Cumberland 
Coast 

English 
Channel 

h 

16  M 

24.5  M 

28.5  M 

T 

44,712  sec 

44,712  sec 

44,712  sec 

"max 

0.59  M/sec 

0.4  M/sec 

- 

^'">T/p 

- 

- 

0.51  M/sec 

<K>J 

0.003  MVsec 

0.0006  MVsec 

- 

\z 

- 

- 

1477  sec 

<h>l 

(observed) 

50  MVsec 

50  MVsec 

18  MVsec 

<^x>T 

(predicted) 

50.4  MVsec 

8.6  MVsec 

9.3  MVsec 
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The  values  of  <E  >j   listed  for  Liverpool  Bay  and  the  Cumberland 
Coast  were  determined  by  Bowden  from  velocity  and  salinity  profiles 
taken  at  one  hour  intervals  throughout  the  tidal  period.  From  these 


he  was  able  to  determine  the  dispersive  mass  transport  <u"  c">y.  This 
was  then  compared  with  the  horizontal  salinity  gradient  through  a  rela- 
tionship such  as  that  given  by  (III-54)  to  yield  a  value  for 
<E  >j.  Bowden 's  reported  values  of  <K  >y  were  determined  by  assuming 
that  the  vertical  distribution  of  salinity,  averaged  over  a  tidal 
period,  was  in  a  steady  state  condition  in  the  oscillatory  sense  due 
to  the  combined  effects  of  horizontal  advection  and  vertical  diffusion. 
Then  using  his  observed  values  of  u  and  c  he  was  able  to  compute 
<K  >j.     Bowden's  results  show  that  in  each  case  <K  >j  varies  approxi- 
mately parabolical ly  with  depth,  being  a  minimum  near  the  surface, 
reaching  a  maximum  near  mid-depth,  and  decreasing  toward  the  bottom. 
It  has  been  assumed  here,  as  previously  stated,  that  this  parabolic 
distribution  may  be  effectively  replaced  by  a  constant  value  corres- 
ponding to  its  mean  over  depth. 

Bowles  et  al.  (30)  made  no  attempt  to  obtain  experimental  values 
for  K  .  However,  they  did  provide  suitable  data  for  the  determination 
of  T  directly.  In  a  series  of  8  experimental  runs  in  which  dye  v;as 
released  as  an  instantaneous  source  at  the  sea  bed  they  experimentally 
determined  the  distance  downstream  at  which  the  dye  was  considered  to  be 
mixed  over  the  water  column.  Thus,  knowing  the  depth  mean  velocity  one 
can  estimate  the  vertical  mixing  time.  Values  of  u„^^  were  obtained  from 
<u>-r   usinn  the  temporal  and  spatial  adjustment  presented  in  (V-5)  on  page 
111. 
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The  comparison  of  results  between  the  observed  and  predicted  values 
for<E>,  is  illustrated  by  Figure  20.  The  observed  value  of  <E^>^  for 
the  Cumberland  Coast  corresponds  to  E'  =  7  x  10"^  which  plots  off-scale 
and  is  indicated  by  an  arrow.  As  shown  in  Table  3  this  value  is  slightly 
less  than  6  times  the  predicted  value.  However,  this  is  not  surprising 
if  one  considers  the  physical  difference  between  the  Cumberland  Coast  and 
an  infinitely  wide  rectangular  channel.  In  fact,  the  agreement  between 
the  predicted  and  observed  values  of  < E  >y  for  the  remaining  two  locations 
is  remarkable.  For  Liverpool  Bay  the  two  values  are  essentially  the  same, 
whereas  for  the  more  open  English  Channel  they  agree  within  a  factor  of  2. 

The  most  comprehensive  set  of  data  available  for  evaluating  values  of 
<E^j   predicted  using  the  methods  presented  here  is  given  by  Bowden 
(31,  5),  and  Fischer  (12)  for  the  Mersey  Estuary  Narrows.  The  Mersey 
Narrows  has  a  width-to-depth  ratio  of  approximately  67  and  a  T'  value 
of  4.25,  Thus,  from  Table  2  it  is  seen  that  this  case  must  be  treated 
as  a  three  dimensional  oscillatory  shear  flow.  In  reporting  his  data 
Bowden  notes  the  significant  effect  of  vertical  density  stratification 
on  longitudinal  dispersion.  Observed  values  of  the  vertical  eddy  diffu- 
sivity  are  one  to  two  orders  of  magnitude  below  the  expected  value  for  a 
comparable  situation  in  a  homogeneous  fluid.  Bowden  attributes  this  to 
an  observed  top  to  bottom  salinity  difference  of  1.3  to  5.0  parts  per 
thousand  and  concludes  that  the  small  values  of  <K  >j  contribute  signifi- 
cantly to  the  large  observed  values  of  <E  >j.  Table  4  summarizes  the 
observed  data  and  gives  the  predicted  value  of  <E  >j  obtained  from 
(IV-43).  The  value  of  K  given  in  Table  4  was  determined  by  using 
Ward's  (23)  experimental  result  that 
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K    =  0.45  u*h  (V-4) 

y 

TABLE  4     OBSERVED  DISPERSION  DATA  AND  PREDICTED 
<E  >^  VALUES   FOR  THE  MERSEY  NARROWS 
X  T 

Parameter  Value 

h  19.5  M 

w  1300  M 

u  1.41  M/sec 

max 

7  44,712  sec 

u^  0.04  M/sec 

<l  >  0.002  M^/sec 

<y-p,  (Eq.   (V-4))  0.35  MVsec 

T  190,100  sec 

cz 

T  1,203,700  sec 

cy 

T'  4.25 

z 

T"  0.16 

c 

Range  of  Observed  <E^>^  161  to  360  MVsec 

<E  >T  using  <u"c">-r  data  170  MVsec 

X  T  ^  T 

<E  >T  using  salt  balance  265  M^/sec 

X  T  ^ 

<E^>^  predicted  by  Eq.   (IV-43)  259  MVsec 

Several   values  of  <E  >-j.  based  upon  observed  data  were  reported  by  Bowden 
all   of  which  fall  within  a  range  of  161  to  360  MVsec.     Bowden  (5)  obtained 
values  of  170  and  265  MVsec  using  measured  <u"  c">-|.  data  as  discussed 
previously  and  a  steady  state  salinity  analysis  respectively.     These  two 
values  for  <E  >^  are  listed  separately  in  Table  4  in  addition  to  the 
observed  range  for  <E  >j. 
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Figure  21  shows  the  observed  range  for  <E^>^  and  Bowden's  two  values 
obtained  from  observed  data  as  they  compare  to  the  predicted  resonance 
curve  for  E'  in  the  Mersey  Narrows.  The  predicted  curve  passes  almost 
exactly  through  the  mid-point  of  the  observed  range  of  E^  values  and  is 
within  3  per  cent  of  Bowden's  observed  E'^  value  using  the  salt  balance 
analysis.  However,  the  predicted  E^  is  52  per  cent  greater  than  Bowden's 
observed  E'  using  <u"c">^.  In  general,  the  agreement  is  considered  to 

be  good. 

All  of  the  longitudinal  dispersion  data  discussed  thus  far  for 
oscillatory  flows  have  been  obtained  in  the  field.  There  are  surprisingly 
few  good  data  obtained  in  the  laboratory  for  open  channel  oscillatory 
flow.  Fukuoka  (32)  has  conducted  an  extensive  set  of  experiments  designed 
to  look  at  this  problem  but  unfortunately  he  made  no  attempt  to  measure 
K  or  K  .  Segall  (19)  also  conducted  a  series  of  laboratory  experiments 
to  study  longitudinal  dispersion  in  oscillatory  open  channel  flow  in 
which  he  experimentally  determined  a  mean  value  for  K^.  A  surmary  of 
Segall 's  data  along  with  the  predicted  value  for  <E^>^  using 
(IV-43)  is  presented  in  Table  5.  To  determine  a  value  for  T^^,  the 
lateral  eddy  diffusivity  was  estimated  by  using  (V-2)  and  (V-4), 
where  u  is  replaced  by  <u>T/  and  it  is  assumed  that 

<J3.   =_i-  u^,,  (V-5) 

T/    3  IT   max 
'2 
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FIGURE  21     COMPARISON  OF  PREDICTED  AMD  OBSERVED  LONGITUDINAL 
DISPERSION  FOR  THE  MERSEY   ESTUARY 


TABLE  5     OBSERVED  DISPERSION  DATA  AS   REPORTED 
BY  SEGALL  AND  CORRESPONDING  PREDICTED 

<E   >,  VALUES 
X  T 
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Parameter 


max 


<Ky>T 


cz 


cy 


<E  >,  observed 
X  T 


<E  >j  predicted  using  (IV-43) 


Value 
0.44  ft 
3.0  ft 
437  sec 
0.1  ft/sec 

0.001  ftVsec 
0.00084  ftVsec 
193.6  sec 
2678.6  sec 
0.07 

0.013  ftVsec 
0.009  ftVsec 


Figure  22  shows  the  E'   resonance  curve  for  T^  =  0.07  corresponding 
to  Segal! 's  data  with  his  measured  value  of  E^  plotted.     It  is  seen  that 
the  predicted  value  for  E^  is  approximately  70  per  cent  of  the  observed 
value.     It  should  be  noted,  however,   that  as  shown  by  Figure  22  the 
predicted  value  of  E'    is  extremely  sensitive  to  T^  for     the  conditions 
used  by  Segal  1.     Thus,  errors  made  in  the  determination  of  h,   K^,  or  T 
could  have  a  significant  effect  on  the  predicted  value  for<E^>y.     Moreover, 
it  appears  that  Segall's  experimental   value  of  0.013  ftVsec  for  <E^>-|- 
was  determined  prior  to  the  establishment  of  a  true  periodic  condition 
for  the  concentration.     As  discussed  in  II-C  this  condition  is  required 
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FIGURE  22  COMPARISON  OF  PREDICTED  AND  OBSERVED  LONGITUDINAL 
DISPERSION  USING  SEGALL'S  DATA  (19) 
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for  the  existence  of  a  stationary  random  dispersive  process.  This  is 
illustrated  by  Segall's  Figure  6.1  in  which  it  is  seen  that  <E^>^  is 
still  decreasing  with  time  during  tidal  cycles  2,  3,  and  4  which  he 
used  to  determine  his  experimental  coefficient.  However,  the  rate  of 
decrease  becomes  small  by  tidal  cycle  4  and  if  one  extrapolates  the 
data  shown  it  is  seen  that  a  more  correct  value  of  <E^>^  from  Segall's 
experiment  would  be  very  close  to  the  predicted  value  of  0.009  ftVsec 

given  in  Table  5. 

The  importance  of  having  good  estimated  values  for  T^^,  T^^,  and 

T'  has  been  demonstrated  by  the  generally  good  agreement  found  between 
c 

observed  and  predicted  <E  >j   values  when  adequate  data  exist  and  the 
correspondingly  poor  agreement  found  between  these  values  when  adequate 
data  do  not  exist.  This  applies  to  both  unidirectional  and  oscillatory 
shear  flow  systems  but  is  of  particular  importance  to  the  latter  due  to 
its  resonant  characteristic. 


CHAPTER  VI 
SUMMARY  AND  CONCLUSIONS 

Analytical  expressions  for  the  cross-sectional  variation  components 
of  velocity  and  concentration,  u"  and  c",  and  the  longitudinal  dispersion 
coefficient,  E.  and  <E  >j,   have  been  developed  using  Taylor's  (1)  defini- 
tion of  dispersion  for  4  cases  of  shear  flows: 

(1)  Unidirectional  flow  in  an  infinitely  wide  rectangular  channel, 

(2)  Oscillatory  flow  in  an  infinitely  wide  rectangular  channel, 

(3)  Unidirectional  flow  in  a  rectangular  channel  of  finite  width, 

(4)  Oscillatory  flow  in  a  rectangular  channel  of  finite  width. 

It  was  shown  in  Chapter  II  that  Taylor's  definition  of  the  longitudinal 
convective  mass  transport  of  a  substance  with  respect  to  the  mean  flow 
is  valid  in  unidirectional  flows  under  the  assumptions  of  steady  uniform 
flow.  In  oscillatory  flows  the  required  assumptions  are  that  the 
velocity  and  concentration  functions  be  temporally  periodic  and  that  the 
flow  be  uniform.  When  these  requirements  are  satisfied  for  their  respec- 
tive flow  conditions  the  time  averaged  convective  mass  transport  of 
substance  with  respect  to  the  mean  flow,  <u"c">^,  is  a  stationary  random 
process  and  is  therefore  Gaussian  in  behavior.  This  permits  its  reformu- 
lation as  a  Fickian  flux. 

Specific  conclusions  regarding  the  predicted  functional  behavior  of 
longitudinal  dispersive  mass  transport  and  the  comparison  of  the  analytical 
solutions  developed  in  Chapters  III  and  IV  with  observed  data  are: 

(1)  Because  of  the  distinctly  different  hydrodynamic  and  kinematic 
characteristics  of  unidirectional  and  oscillatory  flow  systems  the 
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functional  behavior  of  longitudinal  dispersive  mass  transport  in  each 
system  is  unique.  This  fact  is  most  easily  seen  for  two  dimensional 
shear  flows  in  which  the  dispersion  coefficient  for  unidirectional 
flow  varies  directly  with  the  time  of  vertical  mixing  whereas  for 
oscillatory  flow  the  dispersion  coefficient  is  governed  by  a  resonant 
interaction  between  the  period  of  oscillation  and  the  time  of 
vertical  mixing.  It  is  shown  that  as  the  period  of  oscillation 
becomes  very  large  the  dispersion  coefficient  approaches  one-half 
the  value  of  the  coefficient  for  a  comparable  unidirectional  flow. 

(2)  The  description  of  longitudinal  dispersion  in  oscillatory  flow 
through  its  relationship  to  a  "corresponding"  process  in  unidirec- 
tional flow  is  misleading  for  reasons  stated  above.  Although 
technically  correct  values  for  <E^>j  may  be  obtained  in  this  manner, 
great  care  should  be  exercised  in  the  selection  of  expressions  used 
for  <E  >,  and  E.  .  This  choice  is  by  no  means  arbitrary  and  ideally 
should  be  limited  to  companion  solutions  of  the  steady  state  and 
oscillatory  boundary  value  problems  of  flow  and  transport  diffusion 
such  as  those  presented  in  Chapters  III  and  IV. 

(3)  A  solution  for  the  longitudinal  dispersion  coefficient  in  two 
dimensional  oscillatory  shear  flow  is  developed  and  presented  in  the 
form  of  a  single  non-dimensional  resonance  curve  whose  maximum 
occurs  when  the  vertical  mixing  time  is  1.58  times  the  period  of 
oscillation.  This  curve  may  be  used  to  obtain  predicted  values  of 
the  longitudinal  dispersion  coefficient  for  channels  satisfying 

the  two  dimensional  criteria  discussed  in  Chapter  IV  by  knowing: 
(a)  the  period  of  oscillation,  (b)  the  surface  amplitude  of  the 
velocity,  (c)  the  water  depth,  and  (d)  the  depth  mean  value  of 
vertical  eddy  diffusivity. 
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(4)  The  longitudinal  dispersive  mass  transport  in  a  rectangular 
channel  of  finite  width  is  synnetric  about  T^  =  1  for  both  unidirec- 
tional and  oscillatory  flows.  Thus,  a  channel  whose  half-width  was 
n  times  its  depth  would  exhibit  the  same  longitudinal  dispersion  as 
a  channel  whose  depth  was  n  times  its  half-width  provided  that: 

(a)  u    is  the  same  for  both  cases, 
^  '        max 

(b)  The  vertical  and  lateral  mixing  times  for  the  first 
case  were  equal  respectively  to  the  lateral  and  vertical 
mixing  times  for  the  second  case, 

(c)  The  periods  of  oscillation  for  oscillatory  flows  were 
the  same. 

(5)  The  additional  feature  of  a  finite  channel  width  in  oscillatory 
shear  flow  can  either  increase  or  decrease  the  longitudinal  disper- 
sion. The  sense  of  this  effect  is  determined  by  the  relative 
magnitudes  of  the  vertical  mixing  time,  lateral  mixing  time,  and 
period  of  oscillation.  For  T^  =  1  the  peak  of  the  resonance  curve 

for  <E  >T  occurs  at  T'  (or  symmetrically  T')  =  4.5.  As  T'  approaches 
X  T         z  y  ^ 

zero  the  peak  first  increases  in  magnitude  reaching  a  maximum  for 

T'  =  0.25  and  then  decreases  to  the  limiting  value  for  the  two 
c 

dimensional  oscillatory  shear  case.  At  this  point  the  peak  occurs 
at  T'  =  1.58  and  the  resonance  curve  for  <E  >,  is  identical  to  the 

Z  '^  ' 

one  obtained  for  two  dimensional  oscillatory  flow.  Predictive 
values  for  the  longitudinal  dispersion  coefficient  are  obtainable 
from  these  curves  by  knowing:  (a)  channel  width,  (b)  channel 
depth,  (c)  centerline  surface  amplitude  of  velocity,  (d)  period 
of  oscillation,  and  (e)  cross-sectional  mean  values  of  the  verti- 
cal and  lateral  coefficients  of  eddy  diffusivity. 
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(6)  Longitudinal  dispersive  mass  transport  in  finite  width  rectangu- 
lar channels  with  unidirectional  shear  flow  does  not  approach  the 

two  dimensional  solution  in  the  limit  as  one  of  the  channel  dimensions 
becomes  very  large.  In  fact,  E|_  becomes  infinitely  large  as  the 
channel  geometry  is  skewed  in  either  depth  or  width.  This  supports 
Fischer's  theory  regarding  the  dominant  effect  of  lateral  shear 
on  longitudinal  dispersion  in  natural  channels  where  the  width-to- 
depth  ratio  is  much  larger  than  unity. 

(7)  A  comparison  has  been  made  between  predicted  and  observed 
values  of  the  longitudinal  dispersion  coefficient  using  solutions 
developed  in  Chapters  III  and  IV  and  the  field  and  laboratory 
data  of  several  investigators.  These  comparisons  demonstrate 
the  necessity  for  good  estimates  of  the  vertical  and  lateral 
mixing  times  as  part  of  all  longitudinal  dispersion  experiments. 
Where  these  data  are  available  the  agreement  between  the  pre- 
dicted and  observed  coefficients  is  very  good.  However,  where 

values  of  T   and  T   must  be  estimated  the  data  become  very 
cz     cy 

scattered.  This  is  particularly  apparent  for  the  three  dimensional 
unidirectional  shear  flow  case  for  which  nearly  all  T^^  and  T^^ 
values  had  to  be  estimated. 


APPENDIX  A 

DERIVATION  OF  NORMALIZING  EXPRESSIONS 
FOR  TWO  DIMENSIONAL  OSCILLATORY  SHEAR  FLOW  SOLUTION 

A-1  Surface  Amplitude  of  Velocity,  u^^^ 

To  obtain  an  expression  for  u„^^  as  a  function  of  the  pressure  gra- 

maA 

dient  modulus,  K,  the  solution  for  u(z,t)  is  first  evaluated  at  the  free 
surface  and  then  maximized  with  respect  to  time.  The  desired  expression 
for  u(z,t)  is  given  by  Equation  (III-23)  as 

/  4.>  _  iK  n   cosh  B(l+i)(h-z)-i  „iat         /.  -.x 
u(z.t)  =  -  [1  -  cosh  g(1+i)h  ^  ^  ^'^^^ 


where. 


V  2; 


To  begin,  evaluate  (A-1)  at  z  =  h, 

iat 


where. 


and. 


(A-2) 


"(h.t)  =  ^    [1  -  cosh  B(Ui)h  ^  ^'-'^ 


cosh  B(l+i)h  =  J-  (cosh  23h  +  cos  26h)^  e^°''         (A-4) 


-   4--...-^/Sinh  Bh  sin  Bh^  fn_c\ 

""'  -   *^"  ^cosh  Bh  cos  Bh^  ^^  ^^ 
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Next,  substitute  (A-4)  into(A-3)  and  take  the  real  part  of  the 
resulting  expression  to  yield 

n  r  /t.  a.\T    K  .   .  .  /F  K  sin  (at-aa)  /.  /-% 

R^[u(h,t)j  =  -  -  sin  at  +  ^ — r-      (A-6) 

^  ^        a[cosh  26h  +  cos  26h] " 

The  surface  amplitude  of  the  velocity  is  given  by  (A-6)  evaluated  at 

time  t  =  t   ,  i.e.  the  time  at  which  the  maximum  velocity,  u   ,  occurs, 
max  m<i^ 

To  determine  this,  (A-6)  is  differentiated  with  respect  to  time  and  set 
equal  to  zero.  Thus, 

-fZ  Y.   (cos  ot^   cos  az  +  sin  at    sin  aa) 

„     .    ,  max  Ilia  A _  « 

-K  cos  at^^^  + p 0 

^^^  (cosh  2Bh  +  cos  26h) ' 


Solving  for  at^^^. 


at„,^  =  tan"M  _ 

"1^^       1  /F  sin  a 


(cosh  2Bh  +  cos  2Bh)'^  -  a/T  cos  az    I  ,^_-j^ 


Equation  (A-7)  is  then  substituted  into  (A-6)  for  at  to  obtain  the 

desired  expression  for  u^^^.  To  simplify  the  algebra,  let 

max 


F  =  —  (cosh  2Bh  +  cos  2Bh)'^        (A-8) 


Then, 


F  -  cos  a2 

sin  cft^;,v  =  P      (A-5) 

^^^       (F^  -  2F  cosaa  +  1)^ 
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and. 


sin  a2 

cos  at   =  ^  (A-10) 

'"^^   (F2  -  2F  cosa2  +  D' 


Upon  substitution  of  (A-8)  through  (A-10),  Equation  (A-6)  becomes 

u    = j  \   F  cosa2  -  F^  +  F  cosaa  -  cos^aa  -  si 

•"^^   aF(F2  -  2F  cosaz  +  1)" 


n^aa  [ 


or, 


K 


Vx  =  -^(^  -2Fcosa2-Hl)^  (A-ll) 


Now,  from  (A-5) , 


cosa2  =  ^°^^  ^\   "°^  ^^  (A-12) 


Substituting  (A-12)  into  (A-ll)  and  solving  for  K  yields, 

-  a  F  u 
J,  _  mx 

[cosh  23h  +  cos  23h  -  2  cosh  gh  cos  3h  +  1]'^ 

which,  after  some  algebra  and  the  substitution  for  F  using  (A-8)  yields 

as  a  result, 

^  Lin,=>v  (cosh  2eh  +  cos  26h)^ 

K  =  -  —^ (A-13) 

v2  (cosh  3h  -  cos  Bh) 

Equation  (A-13)  is  seen  to  be  identical  to  (III-41)  by  substituting 

2tt 

•j-  for  a  and  applying  (III-36)  which  states 


3h  =  -An 


then, 


z  '       z' 
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VFtt  u  _  (cosh  2  '/irr  +  COS  2  -/ttTI)^ 

K  = 5^^^^ -^^ — r ^—       (A-14) 

T  (cosh  vmtTI  -  cos  i/ttTI) 


This  is  the  expression  used  in  Chapter  III  to  relate  the  pressure  gradient 
modulus,  K,  to  the  surface  amplitude  of  velocity,  u 


max 
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A-2  Surface  Excursion  Length.  L 

The  expression  of  the  pressure  gradient  modulus  as  a  function  of  the 
surface  excursion  length  is  simplified  through  the  use  of  (A-14).  This  is 
seen  by  first  noting  that  the  surface  excursion  length  is  defined  as  the 
product  of  the  time  mean  surface  velocity  during  one-half  of  a  period  of 
oscillation  and  the  time  of  flow,  or 


L  =  <u(h,t)>  ^1     .   2 


(A-15) 


However,  the  surface  velocity  amplitude,  u^^,  is  given  by 


Thus, 


L  =  1  u  (A-17) 

IT  max 


Combining  (A-17)  with  (A-14)  yields  the  desired  relationship  between 

K  and  L, 

V2  TT^  L  (cosh  2  V^  +  cos  2  aATt!)^' 

K  = —r^ — (A-18) 

T^(cosh  virT'  -  cos  VttT^) 


APPENDIX  B 

MATHEMATICAL  DETAILS  FOR  THREE  DIMENSIONAL 
OSCILLATORY  SHEAR  FLOW  DISPERSION 

B-1  Velocity  Solution 

The  governing  partial  differential  equation  for  the  spatial  distri- 
bution of  velocity  over  the  flow  cross  section  is  given  by  Equation  (IV-3) 
as 


i 


^  +  I.. ^-  i  au  =  K  (B-1) 


^  8z^   y  ay^     ' 


for  which  a  solution  of  the  form 


X]      X^  a„    COS  ^^^Jll^  COS  iiJilla  (B.2) 


u  (y,z)  -  ^: ,  ^ ,  a   ^.uo   2h 

^      m  =  0   n  =  0  *""      ^" 


is  assumed. 

Substituting  (B-2)  into  (B-1)  and  performing  the  implied  operations 
yields, 

m  =  0  n  =  0 


00  00 


y        y      fezl^   2h  ^  ^  £y  L  W^^     ""  ^""^    '  (B-3) 


.      cos  i^m+lKz    ^^^  (2j^  ^  ., 
mn  2h  w 
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Next  multiply  (B-3)  by  cos  ^ 5h^    and  integrate  from  -h  to  h, 

noting  that. 


2h 


h  [   0,  )l  ?f  m 

cos  il^j-Qui    cos  12I1L±21IIZ  dz  = 

-h  L  h,    A  =  m 


Equation  (B-3)  is  then  reduced  to  a  single  sum  in  n  and  may  be  written 
after  all  integrations  have  been  carried  out  as 


^\  ^Hr^y  *  -y  l^^^y  *  ia) 


n  =  0 


a   cos  ^^"^^^^y  =  -  f  ^-;r 

^Zn  ^°^        w  (2  il+l)7T 


(B-4) 


To  eliminate  the  summation  in  n  multiply  (B-4)  by  cos  -^ — '-^ 

and  integrate  from  -w/2  to  w/2,  where  now 
sw/2 


cos 


(2k  +  l)7Ty    ^^^  (2  n+l)uy 


cos 


dy 


■w/2 


"  0     ,         n  }«  k 
.  w/2,        n  =  k 


Equation  (B-4)  thus  becomes 


^-^  [i^lpL]^  .  -e^  l^^^y  ^  ia}    a 


8K  (-1)^ 
w(2  JI+1)tt 


ilk 
w/2 

-w/2 


cos 


iLM}ll.  dy  (B-5) 


Carrying  out  the  remaining  integration,  solving  (B-5)  for  a^^j^,  and 
replacing  Z  with  m  and  k  with  n,  the  solution  for  the  Fourier  coefficient 
is 


16  K  (-1) 


m  +  n 


'""^  ^  "  (2  m+1)(2  n+l)7T^   {I     [^^  ''tl^'^y  ^  I     [^^  ^t'^'^'y  ^  ia} 


2h 


w 
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This  is  seen  to  coincide  with  (IV-9).  The  desired  solution  for  the  velocity 
distribution,  u(y,z,t),  is  then  obtained  by  substituting  (B-6)  into  {B-2) 
and  multiplying  the  result  by  e   ,  or 

V^V^  .16  K(-1)"""  cos  11|^  cos  iJ-gl^  e^°^ 

which  completes  the  required  derivation  of  (IV-5). 


(B-7) 
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B-2  Concentration  Solution 

The  governing  equation  for  the  spatial  distribution  of  the  concentra- 
tion variation  over  the  flow  cross  section  is  given  by  (IV-20)  as 

where. 


CO  00 


u;  (y,z)  =  X]  X]  ^»  f"^  ^^^TR^  "=  -^^^T^  -  a'-"-"))  (B-9) 


m  =  0  n  =  0 


^mn  ^'^  ^^^^"  ^^  ^^"^^'  ^"^ 


S^-"'")  =  u^(2m'li(2n+l)  ^^-10) 

Equation  (B-9)  is  then  used  to  substitute  for  u"  in  (B-S)  and  a 
solution  for  c"  of  the  form 

OO         00 

£  =  0  k  =  0 
is  assumed  and  substituted  into  (B-8)  to  yield 

00  00 

X]  X]  <'^z  (r)  '*  %  (^)  ^  ^  ^•°>  »k  "^  T  -^  ^  = 

5,  =  0    k  =  0 

'-  l>  £    £  =-  {-  ^^IW^    cos  iiJlUM  .  g(.,,]     (B.,2) 
m  =  0     n  =  0 

To  solve  for  a'^^,   (B-12)  is  first  multiplied  by  cos  ^^  and  integrated 

from  -w/2  to  w/2.  On  the  left  hand  side  of  (B-12)  the  integration  is 
.w/2  r  0,     k  ?«  j 


^cos^k^dy-  (B-13) 

-w/2  [w/2,    k  =  j 
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which  reduces  the  left  hand  side  to  a  single  sum  in  z  multiplied  by  w/2. 
However,  the  integration  over  y  on  the  right  hand  side  of  (B-12)  does  not 
cause  a  similar  reduction  in  the  number  of  sums  because  of  the  non-ortho- 
gonality of  the  cosine  terms.  This  is  seen  by  defining  the  integral  over 
y  by  Ij.  Then, 

.w/2 
I^  =    \    cos  i^nllM  cos  ^  dy  (B-14) 

.w/2 

w/2  w/2 

cos  (2n  -  2j+l)  ^  dy  +^       cos   (2n  +  2j+l)  ^  dy 
0  0 


\ 


-  w     rcos  (n-j)TT     ,  cos  (n+j)Tr  -, 
"  TT     ^(2n  -  2j+l)  '"{2n  +  2j+l)    ^ 

TT     M2n  -  2j+l)  ^  (2n  +  2j+l)^ 


or. 


,      _    2  W(-1)"   -J        r__2nll____.  /R    ^c^ 

^^ V ^(2n+l)2  -   (2j)2>  (B-15) 

The  g(m,n)  term  on  the  right  hand  side  of  (B-12)  does  not  survive  the  inte- 
gration since 

w/2  pW/2 

g(m,n)  cos  ^  dy  =  g(m,n)  \    cos  ^  dy  =  0 
-w/2  -'-w/2 

Equation  (B-12)  may  then  be  written  as 


Z  =  0 


OS         oo 
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The  remaining  sum  on  the  right  hand  side  of  (B-16)  is  eliminated  in 
a  similar  manner,  i.e.  by  multiplying  (B-16)  by  cos  -^nd  integrating 
from  -h  to  h.  From  (B-13)  the  left  hand  side  integral  is 


•h  C  0,  a  ^  f 


■h  (  h, 


cos  Ij:^  cos  i^H^dz  =  ^  (B-17) 


Let  I2  represent  the  right  hand  side  integral  over  z.  Then, 

rh 
I2  =  \   cos  llJBl^  cos  ^  dz  (B-18) 

J-h 


and  proceeding  in  a  manner  similar  to  that  done  for  Ij, 

T  -  4h  (-1)"''"^  r  (2m+l)    ,  /„  ^Q^ 

^2 V ^(2m+l)2  -  (2f)2>  ^^^^^ 


With  the  aid  of  (B-17)  and  (B-19),  (B-16)  may  now  be  written  as, 

m  =  0  n  =  0 
^mn  i-^f^^'^~h^   ni+l)(2  n+1 ) 


[{2m+l)2-(2j)2][(2  n+l)2-(2f)2] 

The  desired  solution  for  a',  is  then  obtained  by  solving  (B-20)  for 
ai-  and  replacing  f  with  Z   and  j  with  k.  In  addition,  to  maintain  the 
general  nature  of  the  sunmation  indices,  m  is  replaced  with  p  and  n  is 
replaced  with  q.  Thus,  the  solution  for  a^|^  after  substituting  for  a 
using  (B-6)  is 
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00  oo 


a'     =  (^) 


256  K  (-1) 


-{£+k) 


^z  '  h 


^mk  {^v  ^Ki^v 


—  ZZ 

■y   ^  w   '  ^°^     p  =  0     q  =  0 

1 


:B-21) 


[(2  p^T)^-(2£)^][(2  q+l)^-(2k)n  {E^[^^  P;^h^  +  e,  [^^  ^^^  ^^  "M^ 


The  spatial  variation  of  concentration,  c"(y,z,t),  is  then  written  as 

(B-22) 


CO  oo 


•'(^■^•t)  -  zZ  zl 


=>•     ^r.c  ^■fz     ^^o  2kTTy     „iat 
^£k  ^k  ^°^  —    ^°^       w       ^ 


A  =  0     k  =  0 


where. 


'Jlk 


"0.      k  =  i  =  0 
'   h,      k  =  0  or  £  =  0 
J,      k  ?«  0  and  £  ?«  0 


(B-23; 


The  factor,  y  .  ,  must  be  included  to  properly  account  for  the  Fourier 
coefficients  a'    ,  a^.  ,  and  a'^..     This  is  seen  by  first  considering  the 
definition  of  c"(y,z,t).     As  expressed  by  (1-4)  and(II-7),  the  concentra- 
tion, c,  is  decomposed  into  its  cross-sectional  mean  and  variation  com- 
ponents by 

c  =  c  +  c" 


and  thus  by  definition, 

r h    ^w/2 


c"  dy  dz  =  0 


(B-24) 


'-h   J-w/2 


Thus,  the  £  =  0  and  k  =  0  term  in  (B-22)  must  be  zero  to  satisfy  (B-24) 
and  thus  Pqq  =  0.  For  the  case  where  either  £  or  k  is  zero  it  is  seen  from 
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(B-13)  and  (B-17)  that  the  surviving  term  is  twice  as  large  as  the  values 
shown.  Since  a^  varies  as  the  inverse  of  these, terms,  the  values  of  ygj^ 
and  Pj^Q  must  necessarily  equal  0.5. 

The  expression  given  for  c"(y,z,t)  by  (B-22)  is  equivalent  to  that 
given  by  (IV-24).  To  express  c"  as  given  by  (B-22)  in  the  same  manner  as 
given  by  (IV-24)  it  is  necessary  to:  (1)  substitute  for  a^j^  using  (B-21); 
(2)  introduce  the  definitions  of  a   and  the  vertical  and  lateral  "mixing 
times;  and  (3)  transform  all  complex  factors  to  polar  form.  Take,  for 
example,  the  factor  {K J-^)"  +  Li^)^  +   ia  }  and  recall  the  following 


definitions; 


z^  h  y    w 


27r 
a    -- 


\,  -  (w/2)VK^ 


Then, 

(\  (f)'  ^  ^  (^)'  +  ia}  =  f  (2^    i^^Y  +  2^  (k^)^  +  i   >  (B-25) 

and  similarly. 
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Transfomiing  (B-25)  and  (B-26)  to  polar  form,  substituting  into  (B-21) 
and  putting  the  resulting  expression  for  a^.  into  (B-22)  yields 


c"(y.z.t)  '  (|)  ^  ^  Z^  ^  Z.  (B-27) 

H  =  0     k  =  0  p  =  0  q  =  0 

y,^(-l)-(^^^)cos^  cos^e^-^^ 
[(2  p+l)2-(2Ji)2][(2  q+l)^-(2k)2]  r"!   (p,q)  R^l   (£,k) 

which  is  identical  to  (IV-24). 


134 


B-3  Longitudinal  Dispersive  Mass  Transport 

The  purpose  of  the  following  development  is  to  present  the  mathema- 
tical details  involved  in  performing  the  implied  integrations  of  (IV-30) 
repeated  here  as 

<m>^  =  /{  [     RJu"(y.z.t)]  •  Rg[c"(y,z,t)]  dz  dy  V    (B-28) 

\j.w/2  J-h  / 


Using  (IV-TT)  and  (IV-24)  for  u"  and  c"  respectively,  the  integrand 
of  (B-28)  may  then  be  written  as, 

OO  00  OO  00  OO  00 

^(u")  ■  ^(c")  M- 1)-  5]i^  E  Z  Z  Z  Z  Z 

i  =  0     k  =  0  m  =  0  n  =  0  p  =  0  q  =  0 


yoJ-l)"^""^"^  cos  ^  cos  ^  cose,  cose. 


'Jlk^  '  h  w 


(2  m+l)(2  n+l)[(2  p+1  )^-(2£)2][(2  q+l)2-(2k)2]  R,='(m,n)  Ri^(p,q)  R,^  {£,k) 
.  fcos  i2_mll)ll  cos  iLj}fniL    .  g(n,.n)]  (B-29) 


2h 
where. 


9i  =  at  -  ai(m,n) 

02  =  at  -  ai(p,q)  -  aaC^.k) 

It  is  first  noted  that  because  Mqq   =  0  is  required  for  zero  contribution 
to  c"  when  both  a   and  k  are  zero,  then  the  term  g(m,n)  will  not  survive 
the  double  integration  over  y  and  z  in  (B-28)  since 
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pW/2    ph 

\      \   g(m.n)  cos  ^  cos  ^  dy  dz  =  0 

J-vi/9        J-h 


if  both  i  and  k  are  not  simultaneously  zero.  Therefore,  g(m,n)  will  be 
dropped  from  (B-29)  and  will  no  longer  be  considered.  The  integrations 
over  depth  and  width  in  (B-28)  are  then 

JiTTZ  ^„,  (2  m+l)TTZ  .,  _  4h  (-1)"^"^  (2  m+1)     ,r  ^^x 
cos  -[—  COS  -^ — 55—^ —  dz  =  • — ^ — ' ^ [a-iU] 

_^  ^        2^         Tr[{2  m+l)^-(2Jl)^] 

r^'  cosilJliiM  cos  21^  dy  =  ^(-1)"''(^"^1)   (B-31) 

using  the  results  of  the  previous  section.  The  instantaneous  mass  trans- 
port at  any  given  time  during  the  period  of  oscillation  is  obtained  by 
substituting  (B-30)  and  (B-31)  into  (B-29),  or 


.   ,  9cx  .  4096  wh  K^^ 
m  =  (-  ^ 


CO  00  00  OO  00  oo 

^       i!,  =  0k  =  0m  =  0n  =  0p  =  0q  =  0 

cos9i  COS92  Ri{l,k)  (B-32) 

6i(£,k,m,n,p,q) 

It  should  be  noted  that  (B-32)  could  be  used  to  obtain  an  expression  for 
the  dispersion  coefficient  as  a  function  of  time.  To  time  average  (B-32) 
the  remaining  cosine  terms  are  expanded  as  follows: 

cosOi  COS62  =  cos  (at-ai(m,n))cos  (at-ai(p,q)-a2(£,k)) 


136 


or 


■cosa2(5,,k)[cos  (2at-ai(m,n)-ai(p,q))  +  cos  (ai(m,n)-ai(p,q)  )]^ 

{B-33 
+  sina2(i,k)[sin{2at-ai(m,n)-ai(p,q))  +  sin(ai(m,n)-ai(p,q))] j 

The  time  average  over  the  period  is  then, 

Y    I  cosBi  cose2  dt  =1  [  cosa2(il,k)  cos  X  +  sina2(il,k)  sin  x\      (B-34) 

Equation  (B-34)  may  be  simplified  somewhat  by  noting  that 


,k 


cosazCil.k)  =  V     (6-35) 


sinazCil.k)  =  -rr^ {B-36) 

R|  (ii.k) 


Thus,  applying  (B-35)  and  (B-36)  to  (B-34)  and  substituting  the 
result  into  the  time  averaged  form  of  (B-32)  the  expression  for 


<iti>y  is 

9c^   2048  wh  K^T^ 

dE,'                    Till 

00 

oo 

oo 

00 

oo 

00 

<rti>^  =  (- 

>: 

>: 

>: 

>: 

>: 

>: 

1  =  0    k  =  0  m  =  0  n  =  0  p  =  0  q  =  0 

^ik   {(R2(<^,k)-1)^  cos  X   +  sin  X] 

6i(Jl,k,m,n,p,q)  (B-37) 

Equation  (B-37)  is  identical  to  (IV-31). 


APPENDIX  C 

DERIVATION  OF  u  ^^  NORMALIZING  EXPRESSION 
max 

FOR  THREE  DIMENSIONAL  SHEAR  FLOW  SOLUTION 


To  obtain  a  relationship  between  the  pressure  gradient  modulus,  K, 
the  centerline  surface  amplitude 
is  first  expressed  in  polar  form  as 


and  the  centerline  surface  amplitude  of  the  velocity,  u  ^  ,  Equation  (IV-5) 

max 


CO         oo 


m  =  0  n  =  0 


(-1)'""^"  cos  ^^  "i+1)ttz  ^q^  (2  n+l)7Ty  gi(at-ai(m,n) ) 


(2  m+l)(2  n+1)  Rf(m,n) 


Equation  (C-1)  is  then  evaluated  at  y  =  z  =  0  and  the  real  part  taken. 


Re["(0.0,t)]  =  -  ^  ^/_^  ^/_^    a^^  cos(at-ai(m,n))         (C-2) 
m  =  0  n  =  0 


where, 

^mn   ; (C-3) 


(2  m+l)(2  n+1)  Rf(m,n) 


The  velocity  as  given  by  (C-2)  will  equal  u   when  it  is  maximized 

max 


with  respect  to  time.  Thus, 
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138 


CO  00 


^    Re[u(0,0,t)]      =  0  = 


zz 


^'     m  =  0     n  =  0    '"^ 


Expanding  the  composite  sine  term. 


^  /         y>      a^^^  j^sin  at  cosai(m,n)  -  cos  cjt  sinai(m,n)|    =  0 


m  =  0    n  =  0 


and. 


CO  00 


CO  CO 


iin  at  y         /   ^  a^^  cosai(m,n)  =  cos  at  2_ ,    /      .  ^mn  sinai(ni,n) 


m  =  0    n  =  0 


m  =  0     n  =  0 


The  phase  angle  at  which  the  maximum  velocity  occurs  is  then 

OO  00 


m  =  0    n  =  0 


at„^„  =  tan      , 
max  \      °°  °° 


X]    X]  'mn  cosa,(m.n) 


Next  define. 


00  00 


[m  =  0    n  =  0 


Then, 


m  =  0    n  =  0 


a^^^  cosai{m,n) 


00  OO 


(C-4) 


OO  OO 


(_m  =  0     n  =  0 


\n  s^'""!^'"'")! 


(C-5) 


y~^.     ^.a^„sina,(m,n) 


sin  at 


m  =  0     n  =  0 


max 


and. 


00  OO 


X]    ^"mn^^^"^^'"'"^ 


cos  at 


m  =  0    n  =  0 


max 


(C-6) 


(C-7) 
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Equation  (C-2)  evaluated  at  t  =  t        becomes. 


8KT 


max  TT^ 


00  00 


"^  ^W    ZZ   Z^    'mn  coscxi(m.n) 


m  =  0    n  =  0 


+  sin  at. 


...-    Z-j    ZZ    'mn  sincci(m,n) 
m  =  0     n  =  0 

Substituting  (C-6)  and  (C-7)   into  (C-8) 


_      8KT 

max      ~  ii3(i) 


OO  00 

m  =  0     n  =  0 


^mn  cosai(m,n) 


+  { 


y.    X     .    ^mn  sinai(m,n) 


m  =  0    n  =  0 


and. 


max  ir3 


8KT     f^. 

Tr3  \    A    ' 


{C-8) 


:c-9) 


The  desired  relationship  between  K  and  u         is  then 


K  = 


U  TT- 

max 


8T<{. 


(C-10) 
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